ANALYTIC REGULARITY FOR LINEAR ELLIPTIC SYSTEMS IN POLYGONS 

AND POLYHEDRA 



MARTIN COSTABEL, MONIQUE DAUGE AND SERGE NICAISE 



Abstract. We prove weighted anisotropic analytic estimates for solutions of second order 
elliptic boundary value problems in polyhedra. The weighted analytic classes which we 
use are the same as those introduced by Guo in 1993 in view of establishing exponential 
convergence for hp finite element methods in polyhedra. We first give a simple proof of the 
known weighted analytic regularity in a polygon, relying on a new formulation of elliptic 
a priori estimates in smooth domains with analytic control of derivatives. The technique is 
based on dyadic partitions near the corners. This technique can successfully be extended to 
polyhedra, providing isotropic analytic regularity. This is not optimal, because it does not 
take advantage of the full regularity along the edges. We combine it with a nested open set 
technique to obtain the desired three-dimensional anisotropic analytic regularity result. Our 
proofs are global and do not require the analysis of singular functions. 
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Introduction 



Motivation. Solutions of elliptic boundary value problems with analytic data are analytic. 
This classical result has played an important role in the analysis of harmonic functions since 
Cauchy's time and in the analysis of more general elliptic problems since Hilbert formulated 
it as his 19th problem. Hilbert's problem for second order nonlinear problems in variational 
form in two variables was solved by Bernstein in 1904 [9]. After this, many techniques 
were developed for proving analyticity, culminating in the 1957 paper [41] by Morrey and 
Nirenberg on linear problems, where Agmon's elliptic regularity estimates in nested open 
sets were refined to get Cauchy-type analytic estimates, both in the interior of a domain and 
near analytic parts of its boundary. 

Analyticity means exponentially fast approximation by polynomials, and therefore it 
plays an important role in numerical analysis, too. Analytic estimates have gained a renewed 
interest through the development of the p and hp versions of the finite element method by 
Babuska and others. In this context, applications often involve boundaries that are not glob- 
ally analytic, but only piecewise analytic due to the presence of comers and edges, and 
therefore global elliptic regularity results cannot be used directly. 

Elliptic boundary value problems in domains with comers and edges have been inves- 
tigated by many authors. Let us quote the pioneering papers of Kondrat'ev [29] and of 
Maz'ya and Plamenevskii [32, 33, 34, 35]. In these works, the regularity of the solution 
and its singular behavior near edges and corners is described in terms of weighted Sobolev 
spaces. Besides their own theoretical interest, these results are the basis for the convergence 
analysis of finite element approximations of the boundary value problems. 

But whereas these results on elliptic regularity of finite order allow to prove optimal 
convergence estimates for the h version or the p version of the finite element method, they 
are not sufficient for proving the (numerically observed) exponential convergence rate of 
the /ip-version of the finite element method. Indeed, as has been shown for two-dimensional 
problems by Babuska and Guo in [4, 5], the convergence analysis of the hp-YEM requires 
the introduction of weighted spaces with analytic-type control of all derivatives, so-called 
"countably normed spaces". Babuska and Guo proved corresponding weighted analytic 
regularity results for several model problems [4, 5, 22, 25]. 

In three-dimensional domains, as soon as edges are present, there is higher regularity in 
the direction along the edge, and in the /ip-version one introduces anisotropic refinement. 
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performed only in the direction transverse to the edge [6]. The corresponding weighted 
spaces have to take this anisotropy into account. In [23, 24] Babuska and Guo have started 
proving estimates in such spaces in a model situation. 

For three-dimensional polyhedra (containing edges and comers) Guo has introduced the 
corresponding relevant spaces in 1993 [21]: The anisotropy along edges has to be combined 
with the distance to comers in a specific way. Since that time, the proof that the regularity of 
solutions of elliptic boundary value problems with analytic right hand sides is described by 
these spaces has been an open problem, even for the simplest cases of the Laplace equation 
with Dirichlet or Neumann boundary conditions. In the error analysis of hp-FEM, such 
regularity estimates have been taken as an assumption [21, 26, 43]. 

In this paper, we first give a simple proof of the 2D weighted analytic regularity result on 
polygons, for Dirichlet and Neumann conditions, using a dyadic partition technique. Then, 
relying on a nested open set technique, we prove anisotropic regularity along edges in the 
framework of the anisotropic weighted spaces introduced and used in [12, 13], but now with 
analytic-type estimates for all derivatives. Combining the previous two steps with a 3D 
dyadic partition technique at polyhedral corners, we obtain the desired analytic weighted 
regularity in a 3D polyhedron. 

We use two types of weighted spaces of analytic functions. The first type is constmcted 
from weighted Sobolev spaces of Kondrat'ev type. These spaces with "homogeneous norms" 
are suitable for the description of the regularity in the presence of Dirichlet boundary con- 
ditions. For Neumann conditions, a new class of weighted analytic function spaces, con- 
structed from Maz'ya-Plamenevskii-type weighted Sobolev spaces with "non-homogeneous 
norms", has to be used. 

It is important to notice that the above spaces naturally contain the singular parts of so- 
lutions, and give an accurate account of their generic regularity. Thus, in contrast with in- 
vestigations such as [27], we do not need to address separately vertex, edge and edge- vertex 
singularities. Our estimates cover regular and singular parts at the same time. 

Analytic regularity estimates consist of regularity estimates of arbitrary finite order in 
which the dependency of the constants on the order is controlled in a Cauchy-type manner. 
The results of this paper contain therefore, in particular, finite regularity estimates of any 
order in anisotropic weighted Sobolev spaces. For polyhedra, these finite regularity results 
are also new in this generality. In particular, our proof covers the statements formulated in 
[12, sec. 3], and our results generalize those of [2, 8]. 

Our proof of analytic regularity estimates is modular in the sense that it starts from low- 
regularity a-priori estimates on smooth domains and proceeds to singular points, edges, and 
finally polyhedral corners by employing the two techniques of dyadic partitions and nested 
open sets. In order to avoid drowning this clear structure in too many technical difficulties, 
we mainly restrict ourselves to the situation of homogeneous elliptic equations with constant 
coefficients. Generalizations to operators with lower order terms and variable coefficients 
will be briefly indicated. They will be discussed in more detail in our forthcoming book 
[16]. 

Organization of the paper. The main results of the paper are Theorems 6.14 and 7.8 in 
sections 6 and 7. The hypotheses of these theorems as well as the definitions 6.10 and 6.11 
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of the relevant function spaces are necessarily rather complicated. To facilitate their un- 
derstanding, we begin in section 1 by illustrating the general results in the classical simple 
examples of the Dirichlet and Neumann boundary value problems for the Laplace equa- 
tion. We also give simplified definitions for the function spaces in the case of a rectangular 
polyhedron where all edges are parallel to the coordinate axes. 

We start the analysis in section 2 by quoting from [17] an elliptic a priori estimate on 
smooth domains with analytic control of derivatives. This estimate improves the readability 
and efficiency of classical proofs of analytic regularity in smooth domains as can be found 
in [41, 40, 31]. We then prove a refinement of this estimate in view of tackling problems of 
Neumann type. 

In section 3, we make use of a dyadic partition technique to construct weighted analytic 
estimates in plane sectors. This technique has been used in a similar framework in [10] for 
weighted Gevrey regularity. It has been first employed for comer domains in [29], then for 
domains with edges [34], and even for the Laplace operator on a polygon with non-linear 
boundary conditions [28]. The technique of dyadic partitions is a powerful tool to prove 
what we call natural regularity shift results near corners. This expression means that from 
two ingredients, namely basic regularity, i.e. a certain weighted Sobolev regularity of low 
order, of the solution, and improved regularity, i.e. high order weighted Sobolev regularity 
or weighted analytic regularity, of the right hand side, one deduces improved regularity with 
the same weight of the solution. This type of regularity result requires very few hypotheses 
on the weight exponents, none at all in the class of spaces (K and A) with homogeneous 
norms and only a bound from below in the class of spaces (J and B) with nonhomogeneous 
norms. 

In section 4, we combine the local estimates to obtain the natural regularity shift for 
polygons in analytic weighted spaces. 

In section 5 we start the three-dimensional investigation with estimates along an edge. 
The fact that there is additional regularity along the edge has been known and studied for 
a long time (see [19, Theorem 16.13], [36, Satz 3.1], [14, Theorem 4.4]). We therefore 
introduce anisotropic weighted spaces in which derivatives along the direction of the edge 
are less singular than transversal derivatives. There are again two classes of spaces - with 
homogeneous norm (spaces M) and with non homogeneous norm (spaces N). Under the 
assumption of a certain local a priori estimate of low order in the neighborhood of an edge 
point, we prove local analytic anisotropic regularity shift along this edge, by combining 
dyadic partition technique and the classical (and delicate) tool of nested open sets. 

In section 6, we treat polyhedral corners. Relying on suitable definitions of various fam- 
ilies of weighted spaces (similar to [38], but with anisotropy along edges), we are able to 
prove the analytic regularity shift for polyhedra by dyadic partitions around each comer of 
a polyhedron. 

In section 7, we prove the main analytic regularity result for solutions of problems given 
in variational form. The proof combines our analytic regularity shift results with known 
results on basic regularity for the variational solutions. On polygons, we use for this purpose 
Kondrat'ev's classical regularity results in weighted Sobolev spaces, and on polyhedra, we 
use recent regularity results by Maz'ya and Rossmann [38]. In this way, we finally obtain 
the weighted analytic regularity of variational solutions in the right functional classes of 
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[21]. For polygons, we thus prove again in a different and simpler way results which were 
first established by Babuska and Guo [4, 22]. For polyhedra, the results are new. 

We conclude our paper in section 8 by discussing various generalizations. For our proofs, 
we choose in this paper the simplest possible framework of second order homogeneous sys- 
tems with constant coefficients and zero boundary data on domains with piecewise straight 
or plane boundaries. In dimension 2, it is a mere technicality to generalize these proofs to 
the case of second order elliptic systems with analytic coefficients and non-zero boundary 
data. In dimension 3, the possible variation of coefficients along edges introduces more seri- 
ous complications and would require to estimate commutators in a systematic way as in [17, 
Lemmas 1.6.2 & 2.6.2]. In comparison, the generalization to homogeneous transmission 
problems with constant coefficients on a polyhedral partition would be much less difficult. 
Whereas the Stokes system could be considered similarly, things are different for regular- 
ized harmonic Maxwell equations, for which it is necessary to detach the first singularity if 
one wants to obtain a valuable result, see [18] in dimension two. 

We denote by H"^{fl) the usual Hilbert Sobolev space of exponent m, by || ■ \\m;n and 
I ■ |m;f7 its norm and semi-norm. The L^(r2)-norm is denoted by || ■ ||o;n or simply by || ■ ||n • 
Boldface letters like H™(r2) indicate spaces of vector functions. 



In this section, we explain the main definitions and results that culminate in Theorems 
6.14 and 7.8, for a simple class of geometrical configurations and for the simplest elliptic 
boundary value problems, namely the Dirichlet and Neumann problems for the Laplace 
equation. 

1.1. Dirichlet conditions. Let us consider the Dirichlet problem for the Laplace operator 
on a domain Q 



for right hand side / E H~^(n). There exists a unique solution u E H^(f2), owing to the 
Lax-Milgram lemma applied to a variational formulation that is strongly elliptic on HQ(r2). 

If n has a smooth boundary, there holds what can be called the elliptic regularity shift: 
For any natural number n, if f E H" (r2), then u E H"+^(f2). Moreover, if the boundary is 
analytic and / belongs to the class A(r2) of functions analytic up to the boundary of fi, then 
u E A(fi). 

a. Polygons. If has a polygonal boundary, the situation is quite different: If for instance 
has a non-convex angle, the solution u does not belong to in general when / G . 

Instead there hold expansions in regular and singular parts: If / is smooth, then for any 
natural number n we can write 



1. Illustration 



(1.1) 
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in f2, 
on d^l , 



(1.2) 
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Here is the set of the corners c of Let uj^ be the angle of Vl at the corner c. Each corner 
singular part has the form' 



(1-3) Wc,n = Xc(?^c) (ic,fcrc"'sin 



The cut-off function Xc, the polar coordinates (rc, 6*0) and the coefficients (ic,fe are related to 
c. The regularity implication 

(1.4) u G W^iSl) and / G H"(fi) ^ m G H"+2((]) 

holds only if n + 1 < ^ for all c. This precludes any regularity in the analytic class. 

It is known since Kondrat'ev that the use of weighted Sobolev spaces allows a better 
description of the regularity of solutions. We introduce now the Kondrat'ev spaces with a 
notation of our own — which facilitates the definition of weighted analytic classes. With 
rc = rr(x) the distance function to the corner c and ^ = (/3c)cg<r ^ ^^'^ ^ weight multi- 
exponent we define the weighted semi-norms 

(1-5) HK;.,,;.= {E||(n-^'"')^>||o.j'' 

\a\=k c€V 

The space K^(f^) is the space of distributions u such that the norm 

m _ 1 

\\u\' 



{EI^Ik;^,^;^}' infinite- 



fc=0 

These norms are qualified as homogeneous because of the shift +|a| for the weight ex- 
ponent, which makes each term homogeneous with respect to dilations with center in the 
corresponding comer. In the case of homogeneous Dirichlet conditions as in problem (1.1), 
an angular Poincare inequality allows to establish the estimate 

(1-6) ll(n^c-^)"llo;n<^Hl;n' ^eHj((]), 

whence the embedding 

(1.7) Hlin) c KUn). 

This is one of the reasons why the K spaces are appropriate for describing the regularity of 
Dirichlet solutions. 

Kondrat'ev's result for problem (1.1) can be phrased as follows 

Theorem 1.1. [29, section 5.4] If the following condition holds for the polygon fl and the 
exponents (3c 

(1.8) 0<-/3c-l<— VcG^ 

then for any natural number n the solution of problem (1.1) satisfies the regularity result: 

(1.9) u G H^(fi) and f G K|+2(^) =^ K|+^(fi). 



^When ^ is an integer, there is a logarithmic term instead of Tc'' sin ( 
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The analytic class that we associate with the family of semi-norms |n| is 

i\ 5 fC .fj 5 i I 

(1.10) A^(fi) = { n G fl K^(f]) : 3C > 0,Vm G N, < C^+^ml}, 

m>0 

and our regularity result is the following. 

Theorem 1.2 (see Thm 7.1 for the general case). If condition (1.8) holds for the polygon Vt 
and the exponents (3^ then the solution of problem (1.1) satisfies the regularity result: 

(1.11) u^V\^{n) and f ^ kp_+2{n) =^ ueA^{n). 

Our proof consists of the combination of Theorem 1 . 1 with the proof of what we call 
natural regularity shift which holds for any weight exponent without limitation 

Theorem 1.3 (see Thm 4.4 for the general case). For any multi-exponent (3 the following 
regularity result holds for solutions of problem (1.1) in the polygon VL 

(1.12) u e K\{n) and f e kg+2{^) =^ u G A^(fi). 

This theorem is "simply" the analytic version of the well-known regularity shift result 

(1.13) ue K^(fi) and / G K^+2(^) =^ ue K^+^i^) 

valid for any n > and any ^. In fact, the proof of Theorem 1.3 consists in showing that the 
norm estimates corresponding to (1.13) are uniform of Cauchy type in the order n, namely 
there exists a constant C > such that for all integer k >2 and all solutions of (1.1) 

The proof of this family of estimates uses local dyadic partitions and Cauchy type estimates 
for smooth domains with analytic boundary, see Theorem 3.2. 

From this short introduction, we see that the expression "regularity result", which gen- 
erally means the existence of estimates for derivatives of the solution, involves a triple of 
function spaces (Uq, F, U) and states the implication 

(1.15) MGUoand/GF =^ u G U. 

In principle U is in a certain sense optimal with respect to F. Comparing the assumptions 
of Theorems 1.1 and 1.3, we can see the important role of the different choices of the basic 
regularity Uq. In Theorem 1.3, the natural regularity shift result is based on a space Uq = 
K^(f2) that has the same weight exponent § as the space U of improved regularity, which is 
K^+^(fi) in (1.13) and in (1.12). In Theorem 1.1, the regularity result for the solution 

of the variational problem is based on the choice of the energy space ((]) for the space of 
basic regularity Uq. 
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b. Polyhedra. If has a polyhedral boundary, at first glance the situation is similar. Besides 
the set of comers, we have the set of the edges e and the distance functions to each 
edge e. For ^ = {I3c}c£v U {(ie}e&s the weighted semi-norms are defined as 

('■'« |"Ik.m.«={e l|{n''^*'}{n©'"'°'}9>llL}'. 

\a\=k cee eeS ^ 

where r^g denotes the distance function to the set ^ of comers (note that r^g is equivalent 
to the product Hcg't'^c '^^ "^^^ space K™(r2) is the space of distributions u such that 
the sum XI ^=0 I^Ik n finite. The conditions leading to a statement corresponding to 
Theorem 1 . 1 depend on 

i) The opening uo^ of the dihedral angle tangent to ^2 along the edge e, 

ii) The Dirichlet limiting exponent A^"^ at the corner c defined as 



(1-17) A,°'^ = -i + ^^,Di; + i 

where is the first eigenvalue of the Laplace-Beltrami operator with Dirichlet 
conditions on the spherical cap Gc cut out by the cone Fc tangent to f2 at c. 

Theorem 1.4. [37, Thm 2] If the following condition holds for the polyhedron VL and the 
exponents § 

(1.18) 0<-/3e-l<— VeG^ and - | < -/3c - | < A^'^ \/c e "rf 
then for any natural number n the solution of problem (1.1) satisfies the regularity result: 

(1.19) u e and f G K^+2(^) =^ K|+^(f^). 

By some direct extension to the technique leading to Theorem 1.2 we could prove the 
corresponding statement in the analytic classes associated with the seminorms \u\ ^ . 
But such a result would not be of any use for proving exponential convergence of finite 
element approximations, which is one of the main motivations for requiring such analytic 
regularity estimates, cf. [6]. To achieve exponential convergence, one has to use anisotropic 
mesh refinements along the edges, and for the proof one needs corresponding anisotropic 
norm estimates. 

Fortunately, tangential regularity along edges holds, which allows (with a little additional 
effort) to work in anisotropic weighted spaces. To simplify this preliminary exposition, we 
assume that all edges are parallel to one of the Cartesian axes (this is the case if is a cube, 
a thick L- shaped domain or a Fichera comer). Then for each edge e we particularize the 

derivatives in the directions transverse 9x ^ or parallel to that edge so that 

d-^ = d^^ a"" . 

X X X 

We define the anisotropic weighted semi-norms 

(1-20) i«iM.M;n={5: ||{n'c**'}{n(^)'-*"}9>lL}'. 

\a\=k cf^'g eg# ^ 
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Note that (1.16) and (1.20) differ by the exponent /3e + |a| replaced with /3e + For a 
general polyhedral domain, the definition is given in (6.9). The space M™(f2) is the space 

of distributions u such that the sum YTk=i) l^livi fc/s n finite and we denote by A^(f2) the 
corresponding analytic class : 

(1.21) kpJSl) = { w e fl M^(fi) : 3C > 0,Vm G N, < C^^'m^]. 

m>0 

The anisotropic weighted analytic regularity result is the following. 

Theorem 1.5 (see Cor 7.9 for the general case). If condition (1.18) holds for the polyhedron 
Vt and the exponents (3, then the solution of problem (1.1) satisfies the regularity result: 

(1.22) u e H\n) and f E A^+2(f^) =^ uE A^(fi). 

Once again our proof consists of the combination of Theorem 1 .4 with the proof of a 
natural anisotropic regularity shift: 

Theorem 1.6 (see Thm 6.8 for the general case). For any multi-exponent (3 such that 

(1.23) 0<-/3e-l and -(3^-1^— VA; G N, Ve G ^ 

the following regularity result holds for solutions of problem (1.1) m the polyhedron Vt 

(1.24) u E K\{n) and f E ^^+2i^) =^ uE A^(fi). 

Note that condition (1.23) — which in the simple situation discussed here is a sufficient 
condition for Assumption 6.5 to hold — is far less restrictive than (1.18). Nevertheless, we 
need this condition for polyhedra whereas no condition at all was required in the polygonal 
case (Theorem 1.3). 

1.2. Neumann conditions. Let us consider now the Neumann problem for the Laplace 
operator on the domain 



(1.25) 



Au = f in 
dr,u = on d^l . 



for an L^(i7) right hand side with zero mean value. There exists a solution u E H^(f2), 
unique up to the addition of a constant. In a smooth domain, u satisfies the elliptic regularity 
shift. 

a. Polygons. The solutions have singular expansions like (1.2) with cos functions instead 
of sin. In particular, independent constants are present at each comer. The scale of weighted 
spaces K^ (f2) cannot yield optimal regularity results because solutions of problem (1.25), 
even regular, do not belong to K^(0) for (3c = —I in general. To overcome this difficulty, 
we consider weighted spaces Ja(fi) with non-homogeneous norms defined as 



(1-26) NL.(.)={E||(n-c^"")^>C} 
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instead of 



NiK;,«,={E||(n''?*"")a?HL}^ 



for K^(f2): the exponent in J spaces does not depend on the order of derivation. As a 
particular case we obtain the standard Sobolev spaces 

H"(fi) = rs^). 

The advantage of this notation is a natural notion of the corresponding analytic class. For a 
multi-exponent (3 we first set 



(1.27) 



Kn := max —/3c 



An important property of the spaces is that if m > kr, then the norm ||n|| is 

equivalent to the following "step-weighted" norm 



(1.28) { 5Z II ( n C'^'^^+l"!'"^) d'^u ^ ^} 



|a|<n. ce'^ 



This implies that we have the continuous embedding of J"^^(r2) into J^(f2) when m > k^, 
which leads to the definition of the analytic class 

(1.29) Bp_{n) = [ue fl f^{n) : 3C>0,Vm>«:^ \M (n) ^ '^■} ■ 

If /3c ^ near the corner c the analytic classes A^(r2) and B^(fi) differ by polynomial 
functions of degree < [—/3c — 1] in Cartesian variables: Such functions are present in B^(r2) 
but not in A^(f2) , see Remark 4.3 . The counterparts of Theorems 1.1 to 1.3 are" 

Theorem 1.7. [30, Thm 7.2.4 and section 7.3.4] If condition (1.8) holds, then for all integer 
n such thatn + 2 > solutions of problem (1.25) satisfy: 

(1.30) u G H^(fi) and f G J|+2(f^) =^ ue J|+^(fi). 

Theorem 1.8 (see Thm 7.3 for the general case). If condition (1.8) holds, then solutions of 
problem (1.25) satisfy: 

(1.31) u G H\n) and f G B^+2(fi) =^ m G B^(fi). 

Theorem 1.9 (see Thm 4.5 for the general case). For any multi-exponent § and any integer 
m > solutions of problem (1.25) satisfy: 

(1.32) u G J^{n) and f G B^+2(^^) =^ ue B^(fi). 



^The fact that the same condition (1.8) works for both Neumann and Dirichlet boundary conditions is a 
very particular case (Laplace operator in 2D). 
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b. Polyhedra. We have to introduce the isotropic and anisotropic weighted spaces with non- 
homogeneous norms J^(f2) and N^(f2) corresponding to K^(fi) and M^(f2), respectively. 
They are defined by their norms as follows: 



ii"ii.j,«={ E ||{^''^"}{^(|)*""K»ILJ^ 

\a\<n c&'g eG<r ^ 

and, when the edges are parallel to the axes: 

(i-st) = { E ||{ n'-r"'-*""^°'}{ n (^)"""'""'-''°'}a?Hi' } 



|al<n cG%f eG# 



for any n> Kg where 



(1.35) 



Kg 



:= max < max —/3c, max — /3e ^ 



Note that the definitions (1.33) and (1.34) are coherent since when n > Kg, the norm of 
J^(fi) is equivalent to the norm obtained by replacing /3c + n by max{/3c + |a|,0} and 
/3e + n by max{/3e + |«|, 0}. 

The anisotropic weighted analytic classes are then defined as 

(1.36) Bgin) = [ue fl N™(fi) : 3C>0,Vm>/t^ H^^H N-(n) ^ ^"^'^'}- 

On the same model as (1.17) we define the 3D Neumann limiting exponent A^^" at the 
comer c as 



(1.37) Aj"" = -- + ,/^"? 



where /i^g" is the second eigenvalue'' of the Laplace-Beltrami operator with Neumann con- 
ditions on the spherical cap Gc- With this we can state the counterparts of Theorem 1.4-1.6. 

Theorem 1.10 (cf [38, Thm 7.1] and [19, Chap 8]). If the following condition holds for the 
polyhedron Q and the exponents ^ 

(1.38) < -/3e - 1< — Ve e ^ and - ^ < -/3c - f < min{2, A^'^} Vc G ^ 

then for any natural number n such that n + 2 > Kg solutions of problem (1.25) satisfy: 

(1.39) u G H^(fi) and f G J|+2(f^) =^ ue J|+^(fi). 

The minimum of Ac with 2 in (1.38) comes from the conditions of injectivity modulo 
polynomials of [19] which replaces the usual spectral conditions when polynomial right 
hand sides are involved, see section 7.2 b. Then the anisotropic weighted analytic regularity 
result in non-homogeneous norms is the following. 



^The first Neumann eigenvalue /^t^|^ is zero. 
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Theorem 1.11 (see Thm 7.1 1 for the general case). If condition (1.38) holds for the polyhe- 
dron and the exponents (3 then solutions of problem (1.25) satisfy: 

(1.40) u e and f G B^+2(f^) =^ ue B^(l]). 

The corresponding natural regularity shift result is the following. 

Theorem 1.12 (see Thm 6.14 and Remark 6.15 for the general case). For any multi-exponent 
(3 such that condition (1.23) holds, we have the following regularity result for solutions of 
problem (1.25) in the polyhedron Vt 

(1.41) (n G J^(fi) with m > k^) and (/ G B^+2(l^)) =^ ue B^(fi). 

Remark 1.13. In Theorem 1.12, condition (1.23) plays the role of a necessary and sufficient 
condition for the validity of Assumption 6.13. The fact that this is the same condition as 
the one in Theorem 1.6 is again a particularity of the Laplace operator in 2D, like we have 
seen for condition (1.8). More precisely, whereas in the Neumann case the lower bound 
< — /3e — 1 is necessary, in the Dirichlet case of Theorem 1.6, the sharp lower bound 
would be — — < —(3e — 1, if we allow stepping out of the variational space H^(r2), see 
§7.2c. 

Remark 1.14. Theorem 1.12 and, if we replace A^^" by A^"^, also Theorems 1.10 and 1.11, 
are true for the solution of the Dirichlet problem (1.1). 

Example 1.15. (i) For the unit cube □ = (0, 1)^, 

— = 2, A°'^ = 3 and A,^^" = 2. 

We 

(ii) For the Fichera comer F = (—1, 1)^ \ (0, 1)^, at its non-convex corner Cq = (0, 0, 0) and 
its non-convex edges e^, z = 1, 2, 3 we have 

— = -, ~ 0.45418 and A,^"' ~ 0.84001. 

We, 3 

Here the 3D Dirichlet and Neumann limiting exponents X^" and A^^^" have been computed 
by Th. Apel using the method of [1]; for the Dirichlet exponent see also [42]. 

1.3. BibUographical comments. On polyhedra, the first proof of a Fredholm theory in type 
K weighted spaces is due to Maz'ya-Plamenevskii [33]. Optimal regularity and Fredholm 
results for coercive variational problems are then proved in [19] using unweighted Sobolev 
spaces H"^ = J™^, in [37] using type K weighted spaces (cf Theorem 1.4), and in [38] 
using type J weighted spaces (cf Theorem 1.10). Let us mention that regularity results in K 
weighted spaces are also proved in [7, 39] by an approach distinct from [37]. 

Concerning finite anisotropic regularity, the first partial result is due to Apel-Nicaise [2] 
for the Laplace-Dirichlet problem; more complete and general statements using spaces M 
are announced in [12] for the same problem; a proof of anisotropic regularity using dis- 
tinct, but similar spaces, is provided in [8]. Finally, concerning analytic weighted regularity, 
prior to the present work, we only find two papers relating to three-dimensional domains: 
[21] where Guo describes the suitable weighted analytic classes of type B, and [24] where 
estimates along edges are given for the Laplace-Dirichlet problem. 
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Recently higher order regularity for the Dirichlet problem on hypercubes was proved in 
a related, but different class of anisotropic weighted Sobolev spaces [20]. This was also 
motivated by questions of convergence of numerical approximations, but with a different 
aim, namely optimal approximations intended to overcome the curse of dimensionality. 



2. Local Cauchy-type estimates in smooth domains 

The starting and key point is a local Cauchy-type ("analytic") estimate in smooth domains 
that is proved by using nested open sets on model problems, like in the Morrey-Nirenberg 
proof [41] of analytic regularity, and a Faa di Bruno formula for local coordinate transfor- 
mations, see [17, Theorem 2.7.1]. This proof, which is given in detail in [17], will not be 
repeated here. 

Proposition 2.1. Let Q be a bounded domain in W^,n>2. Let T be an analytic part of the 
boundary ofQ. Let L be a N x N elliptic system of second order operators with analytic 
coefficients over QUT. Let {T, D} be a set of boundary operators on T of order 1 and 0, 
respectively, with analytic coefficients, satisfying the Shapiro-Lopatinskii covering condition 
with respect to L on T. Let two bounded subdomains Q = U (1 Q and Q' = W Ci Q be given 
with lA and W open in M" and hi C W. We assume that V := d^l' fl d^l is contained in T. 
Then there exists a constant A such that any u G H^(r2) satisfies for all k E N, k > 2, the 
improved a priori estimates ('finite analytic estimates") 



We will know prove a refinement of this estimate where in the right-hand side of (2.1) 
the H^-norm is replaced by the H^-semi-norm. This will be convenient for boundary value 
problems of Neumann type. When L, T and D are homogeneous with constant coefficients, 
this version is a consequence of the previous result, obtained by a simple argument based on 
the Bramble-Hilbert lemma. 

Corollary 2.2. We assume that the operators L, T and D are homogeneous with constant 
coefficients. Let m > 1. There exists a constant A independent of k such that there hold 
the following estimates for all k > m and all u e H'^iVt') satisfying the zero boundary 
conditions Tu = and Du = onT: 

(2.2) —lul ^<A''+^l V -|Lu| +|u| - }. 

l=m.-l 

Proof. We start with any u G and use estimate (2.1). We split the right hand side of 

the inequality into two pieces according to: 



k-2 ^ 
1=0 



(|Lu|,^^,+ ||Tu||^^,p,+ pu||^_^3^p, ) + ||u||^^^, = B*{u) + BM 
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with 

k-2 



i=m—l j=m—l j='m 

£=0 

k-2 ^ min{^,m-2} min{£+l,m— 1} 

Since the orders of L, T and D are 2, 1 and respectively, we obtain 

Since, moreover, the operators L, T and D are homogeneous, we have the invariance of 
B*{u) by subtraction of polynomials of degree less than m — 1 

Altogether, using (2.1) for u — if we obtain for all /c > m 
With the Bramble-Hilbert lemma [11], this gives 

Applying this to functions u satisfying zero boundary conditions, we obtain (2.2). □ 



3. Weighted Cauchy-type estimates in plane sectors 

The model singular domains in two dimensions are the infinite plane sectors. Let JC be 
an infinite sector with vertex at the coordinate origin = (0, 0). In polar coordinates {r,9) 
such a sector can be described as 

(3.1) /C = {x G : uji<9 < W2}, 

where a;2 = 1^1 + w with cui E (— tt, tt), and cu E (0, 2n] is the opening of the sector JC. For 
i = 1,2, let Fj be the side 9 = cuiof JC. 

We consider an elliptic system L in /C and on each side Fj a set of boundary operators 
{Ti, Di} satisfying the covering condition. We assume that the operators L, Ti and Di 
are homogeneous of orders 2, 1 and 0, respectively, with constant coefficients. For any 
subdomain W of JC, we consider the system of local interior and boundary equations 

r Lu = f in /cn W, 

(3.2) i T, u = on F^ n W, i = 1, 2, 

i Au = on FinW, i = 1,2, 

which is the localization to W of the elliptic boundary value problem Lu = f in JC, with 
zero boundary conditions on Fi and F2. 
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3.1. Weighted spaces with homogeneous norms. These spaces coincide with those intro- 
duced by Kondrat'ev in his pioneering study of corner problems [29]. The weight depends 
on the order of the derivatives. We adopt a different convention than [29] in our notation 
in order to make the definition of corresponding analytic classes more natural (see (4.6) 
below). 

Definition 3.1. Let (3 he a real number called the weight exponent, and let m > be an 
integer called the Sobolev exponent. Let W be a subdomain of /C. 

The weighted space with homogeneous norm K™(>V) is defined, with the distance 
r = |x| to the vertex 0, by 

(3.3) K;^(>V) = {ue Lf„,(>V) : r^+l"l9> G Va, |a| < m] 
and endowed with semi-norm and norm respectively defined as 

(3.4) Inl,^ „ = lk^^'"'f^x^*IL i^, and ||n|L„,,^„ = kL , a 

^ ' ' 'K;m,^;W /-^ " * "0;W " "K^(W) A-i' ' 'K;fc,^;>V 

a|=m fc=0 

In these spaces we have the following estimates. 

Theorem 3.2. Let W and W be the intersections ofK. with the balls centered at of radii 
1 and 1 + 6, respectively. Let /9 G M and n G N. Let u G Hfo^(W" \ {0}) be a solution of 
problem (3.2). Then the following implication holds 

(3.5) u G Kj(>V') and f G K;^+2(>^') =^ " e K;+'(W) 

and there exists a constant C > 1 independent of u and n such that for any integer k, 
< k < n + 2, we have 

,3.6) i ( j: \\r^-'"^a:u\\l„f < c-'{ ( E \\r^-^-'''s:nl„,f 

\a\=k i=o \a\=e 

^ II' ^x"ll0;W' / 

|a|<l 

Proof. Let us assume that u G K^(>V') and Lu = f G K^_^^20^')- Let us prove estimate 

(3.6) . By definition of the weighted spaces, the right-hand side of (3.6) is bounded. The 
proof of the estimate is based on a locally finite dyadic covering of W and W". Let us 
introduce the reference annuli, see Fig. 1 

(3.7) V = {x G /C : i < r(x) < 1} and V' = {x e IC : \ - 6 < r(x) <l + 6}. 
and for G N the scaled annuli: 

V^ = 2-^V and V^ = 2'^V. 

We check immediately that 

W = U and W' = U V; . 

Step 1 . We are going to apply Proposition 2.1 in two regions which separate the two sides 
Fi and F2 of /C where the boundary conditions can be distinct. We recall that the sector /C 
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Figure 1 . Reference and scaled annuli for a sector /C of opening 37r/2 

is defined by the angular inequalities ui < 6 < U2. Let 0J3 := ^{ui + ^2)- We define the 
sectors /Ci and IC2 by 

/Ci = {x G : uji<e < Us} and /C2 = {x G : U3 < 6 < U2}. 

Let < 5o < |('^2 — wi). We define the larger sectors IC[ and /C2 by 

/C; = {x G : < < ^3 + 5o} and /C2 = {x G : W3 - < ^ < ^2}- 

Letz G {1,2}. Since the system L is elliptic and covered by its boundary conditions {Ti,Z}j} 
on Fj, the reference domains V fl /C, and V fl /C- satisfy the assumptions of Proposition 2. 1, 
and there exists a positive constant Ai such that for all G N, /c > 2, we have: 

1 ^ 1 

e=o ' £=0 

for any function u satisfying the boundary conditions of (3.2) and f := Lu. From these 
estimates for i = 1, 2 we deduce immediately, with A = max{Ai, A2} 

1=0 ' e=o 

Step 2. Since r(x) ~ 1 on V, we can insert weights in the reference estimate (3.9): There 
exists a positive constant B such that for all A; G N, A; > 2 

(3.10) ^ (E ll''(«)'"""«0;v)' £ ll'-'*)'"'"'"'^?'"";-)' 

|a|=fc £=0 |q|=^ 

+ 5^ ||r(x)^+H5^ti||^^„^|. 

|a|<l 

By the change of variables x — t- x = 2^^x that maps V onto (resp. V onto V/J coupled 
with the change of functions 

u(x) := u(x) and f(x) := Lu which implies f(x) = 2~^''f(x), 
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we deduce from estimate (3.10) that 



\a\=k 

Multiplying this identity by 2"^^"*"^, the above estimate is equivalent to 

1 1 ^'-2 ^ 

-L / X ^ 11 / X 11 2 \ 2 , / X ^ i 



k\ 



\a\=k i=0 \a\=e 

|a|<l 

Summing up the square of this estimate over all /i and considering that only a finite number 
of the V/1 overlap, we get the desired estimate (3.6). □ 

3.2. Weighted spaces with non- homogeneous norms. In these spaces the weight expo- 
nent does not depend on the order of derivatives. Standard unweighted Sobolev spaces are a 
special case. The weighted Sobolev spaces with nonhomogeneous norms allow an accurate 
description of the regularity of functions with non-trivial Taylor expansion at the comers. In 
particular, they are useful for studying variational problems of Neumann type, because the 
variational space does not fit properly into the scale K^. 

Definition 3.3. Let (3 he a real number and m > an integer. 

Let W be an open subset of /C. The weighted space with non-homogeneous norm 
j;^(>V) is defined by 

(3.11) ^J(W) = {ue : r^+'^dy e Va, |a| < m} 

with its norm 

||^||2 ^ \\r^+"'d''uf 

ll"llj™{W) II' "llo;W ■ 

\a\<m 

Note that the semi-norm of J^(W) coincides with the semi-norm of K™(>V). They are 
both denoted by I ■ I „ „ With this notation, we have 

m 

ll"llj!5^(W) ^ l"lK;fc,/3+m-fc; W 

k=0 

Remark 3.4. If W is a finite sector with vertex at the origin, there holds [35, 30] (more 
details are given in [15] and [16, Ch. 11]): 

If /3 > -1, then J^(>V) = K^(>V) for all m e N. 

If /3 < -1 and m < -/3 - 1, then, again, J;^(>V) = K;^(>V). 

If /3 < — 1 and m > —(3 — 1, then one has to distinguish two cases: 
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• the generic case —/3 ^ N, in which one has 

J™(>V) = KJ{W) © P[~''"^] 

where pt^'^^^l is the space of polynomials of degree not exceeding —(3 — 1; 

• the critical case -f3 e N, in which JJ^(>V) contains KJ^(>V) © p-/^-i as a strict 
subspace. 

We deduce from [15, Prop. 3.18] the following "step-weighted" characterization of in 
the case of two space dimensions: 

Proposition 3.5. Let /3 G M and m > —(5 be a natural number. IfW is bounded, then the 
norm in the space J™(>V) is equivalent to 

(3.13) [Y. lk'^^^^'^'"''°^^x«llo;w)'- 

Corollary 3.6. Let /3 eR. Letm> -/3 be a natural number^ Then J™+^(>V) C J™(>V). 

Theorem 3.7. Let W and W be the intersections ofJC with the balls centered at of radii 
1 and 1 + 6, respectively. Let /3 G M. Let m > 1 be an integer such that m + 1 > — /3. Let 
n > m — 1 be another integer Let u G HfQj.(>V' \ {0}) be a solution of problem (3.2). Then 
the following implication holds 

(3.14) u G J^(W) and f G J^+alW) =^ u G J^+'(W) 

and there exists a constant C > 1 independent of u and n such that for all integer k, 
m < k < n + 2, we have 

,3,5) g l{Ey''''"'^<y.'f 

\a\=k e=m-l \a\=e 

\a\=m 

Proof Let n > m - 1 and assume that u G JJ^(W') is such that f G JJ^+2(^')- ^^k = m, 
estimate (3.15) is trivial. So, let k be such that m+l<A;<n + 2. Let us prove estimate 

(3.15) . Since m + 1 > — /3, we have 2+ |a| > —(3 for all a with length > ??7 — 1. Therefore, 
as a consequence of Proposition 3.5, the right-hand side of (3.15) is bounded. 

Then, in a similar way as in the proof of Theorem 3.2, we start from estimate (2.2) written 
for the reference domains V and V' and we apply the same dyadic covering technique. We 
arrive directly at the estimate (3.15). 

It remains to prove that u G J^'''^(W). Since W is bounded, estimate (3.15) implies that 
^/3+n+2^a^, bgiongs to for all a,m<\a\<n + 2. Since u G J^(W), we deduce 

that r'^+^+^S^u also belongs to L^(W) when \a\ < m, which ends the proof. □ 

^ For the sake of simplicity we did not quote [15, Prop. 3.18] in its full optimality. In fact the embedding 
jm+i(-yy) ^ J'»(W) holds as soou as m > - 1. 
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4. Natural weighted regularity shift in polygons 

Let be a polygonal domain. This means that the boundary of is the union of a finite 
number of line segments (the sides T^, for indices s G y). We do not assume that is a 
Lipschitz domain, that is we include the presence of cracks in our analysis. The vertices c 
are the ends of the edges. Let us denote by ^ the set of vertices and 

(4.1) rc(x) = dist(x, c). 

There exists e > such that, setting 

(4.2a) l^c = {x G : Tc < e}, 

we have 

(4.2b) He n He' = 0, Vc ^ c'. 

Setting fic^-* = {x G : Tc < |}, we define 

(4.2c) f]o = f^\|Jf^. 

We also define larger neighborhoods choosing e' > e such that 

(4.2d) n'^ = {xen : rc < e'}, n H^, = 0, Vc 7^ c' , 

and we finally set 

(4.2e) fi[, = fi\|J?f), 

where Q.'^c^ = {x G ^2 : Tc < |}. For each comer there is a plane sector /Cc with vertex 
such that the translation x i-t- x — c sends fic onto /Cc n B{0, e). 

Let ^ = {Pc)cev ^ be a weight multi-exponent and m G N a Sobolev exponent. By 
localization we define the weighted semi-norm on any domain V C fi: 

(4-3) \uf , „ = Y (\\d>\\l ,,r.n + Y lkc''''"'^x W|ln .,r,o 

and the norms, cf. (3.4) and (3.12) 

m m 

(4-4) II"IIk™{V) ^ l^lK;fc,^;V ^'^^ Il^llj™(f7) ^ l^lK;fc,^+m-fc;V ' 

fc=0 ~ k=0 

defining the spaces K^(V) and J^(V), respectively. If all weight exponents /3c coincide with 
the same number (3, these spaces are simply denoted by KJ^'(V) and J™(V), respectively. 
Boldface notations K™(V) and J™(V) indicate vector- valued functions. 

Remark 4.1. The semi-norm 1^1^^ .^^ is equivalent to the globally defined semi-norm 



(4.5) { i|(n'-?"'°')9>||L} 



\a\=m c&V 
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We define on any domain V C ^ tlie corresponding weighted analytic classes. 

(i) With homogeneous norm: 

(4.6) A^(V) = { G n ^I'(^) ■ 3^ > 0'^^ ^ ^' l^l,n,^;v ^ ^'""''^!}- 

■m>0 

( //) With non-homogeneous norm: For a multi-exponent ^ let 

(4.7) K^:=max— /3c. 

As a consequence of Proposition 3.5, for all m > we have the continuous embedding of 
J|*+^ (V) into (V) . We introduce 

(4.8) B^(V) = { n G fl JJ(V) : 3C > 0,Vm > < C^+^ml}. 

Remark 4.2. fz) The classes A^(r2) and B^(r2) can be equivalently defined replacing semi- 
norms |n| o.ci^y global semi-norms (4.5). 

(ii) The classes A^(f2) can also be equivalently defined locally i.e. 

Al3_in) = {ue Lf„,(fi) : e A(fio) and u\^^ G A^^{n,) Vc G ^^}. 

Here A{Qq) is the unweighted class of analytic functions on ^Iq. The spaces B^(f2) allow 
analogous local descriptions. 

Remark 4.3. (i) Our spaces B^(r2) coincide with the family of countably normed spaces 
B^{il), introduced by Babuska and Guo [4]: The spaces B^^(Q) are defined for i E N and 
< /3 < 1, and there holds 

(4.9) Blin) = B^^,{n) . 

(ii) The relation between the classes A^(f2) and B^(fi) follows from the relation between 
the weighted spaces with homogeneous and nonhomogeneous norms K"(f2) and J™(f2). 
As a consequence of Remark 3.4 it follows that for (3 > —1 there holds B/3(r2c) = Ajs^Qc), 
whereas for /3 < — 1 one has in the non-critical case —(3 ^ N: 

(4.10) B^(nc) = A^(l]c)©p[-^-'l 

and in the critical case —(3 G N: B/3(r2c) contains A/3(r2c) © P^^^^ as a strict subspace. 

We consider a "mixed" boundary value problem on the polygonal domain We suppose 
that we are given an homogeneous second order elliptic system L with constant coefficients 
and for each side s a covering set of boundary operators {Tg, D^} of order 1 and 0, homoge- 
neous with constant coefficients. The boundary value problem under consideration is: 

{Lu = f in r2, 
Ts u = on Ts, s G ^ , 
u = on Ts, s G ^ . 

Note that one of Tg or Dg may be the zero operator, in which case the corresponding bound- 
ary condition is empty. 
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We can now prove the following statements of natural regularity shift in weighted analytic 
spaces with homogeneous or non-homogeneous semi-norms, respectively: 

Theorem 4.4. Let (3 = (/3c)c6%f ^ weight multi-exponent. Let u G HfQj,(r2 \^) be a 
solution of problem (4.1 1). For all integer k > 1, there holds the global uniform estimate 

1 ''''^ 1 

1=0 

The following implications hold 

(4.13a) u G K^(fi) and f G K|^2(^) ^ " ^ K^+^Q) (n G N). 

(4.13b) u G K^(n) an J f G A^+2(^^) ^ u G A^(fi). 

Proof. The uniform estimate (3.6) is valid between fic and fi^ for all c G The uniform 
estimate (2.1) of the smooth case is valid between f2o and H'q. Combining these estimates 
we obtain the global uniform estimate (4.12). The finite regularity shift (4.13a) is then 
straightforward. If f G A^+2(r2), it satisfies 1^1^^.^^ < F^~^^i\ for some constant F > 1. 
Thus estimate (4.12) yields 

|ul ,,^<k\ C'+'( y F^+i + ||u|l ,,,, ) = k\ C'^+'f^^^ + ||u|l ,,,, V 

Hence u G A^(f2), which proves (4.13b). □ 

Theorem 4.5. Let (3 = (/3c)ce<g' ^ weight multi-exponent. Let u G Hf^cl^ \^) be a 
solution of problem (4.1 1). Let m > 1 be an integer such that m > —(3c for all c G ^. For 
all integer k > m, there holds the global uniform estimate 

l=m-l 

The following implications hold 

(4.15a) u G J'^in) and f G J^+sl^) ^ " ^ J^+^(l^) (n + 2 > m). 

and 

(4.15b) u G JJ(fi) and f G B^+2(f^) ^ u G B^(fi). 

Proof. The finite regularity shift (4.15a) is an obvious consequence of Theorem 3.7. The 
proof of (4.15b) is similar to that of (4.13b), based on estimate (3.15). □ 
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5. Local anisotropic Cauchy-type estimates in dihedral domains 

Infinite dihedral domains (or wedges) are the model domains for polyhedra which have 
the lowest level of complexity. In this section, we consider dihedral domains V ma model 
configuration, that is there exists a plane sector /C with vertex so that 

(5.1) V = lCxR and x = (xi, X2, X3) = (x^, X3) G P x_l G /C, X3 G M. 

The edge e of the dihedral domain V is the line xi = X2 = 0. 

Let V be any subdomain of V. We consider the system of local interior and boundary 
equations 

{Lu = f in VnV, 
TiU = on(riXM)nV, i = l,2, 
Au = on(riXM)nV, 2 = 1,2, 

where the operators L, Ti and Di are homogeneous with constant coefficients and form 
an elliptic system. The system (5.2) is the localization to V of the elliptic boundary value 
problem Lu = f in V, with zero boundary conditions on Fi x M and r2 x R. 



5.1. Isotropic estimates: natural regularity shift. The weighted spaces for the dihedron 
are defined by the same formulas as in the case of a plane sector: 

Definition 5.1. Let /3 be a real number and let m > be an integer. Let W CV. 

The isotropic weighted spaces K™(>V) and J™(W) are defined, with the distance 
f '■= = i/xf + X2 to the edge e, by 

KJ(W) = {u e iLiW) : r^+^'^^dy e l\W), Va, |a| < m} 
J-(>V) = {ue : r^+'^dy e L2(>V), Va, |a| < m} 

endowed with their natural semi-norms and norms. Recall that denotes the de- 
rivative with respect to the three variables xi, X2, X3. 

We call these spaces isotropic, in opposition with the anisotropic spaces M™(>V) and 
N"(W) which will be introduced in the next subsection. 
We gather in one statement the results concerning the K and the J spaces. Here we set 

w=(/cns(o,i))x(-i,i) 

>V^ = (/Cni3(0, l + e)) X (-1 -e,l + e), e > 0. 

Theorem 5.2. Let /3 eR and n e K Let u e Hf^^(W£ \e) be a solution of problem (5.2) 
with V = We. 



23 



(i) Ifu E K^(We) and f G K|^_^_2(We) u G K^+^(W) anJ ?/zere exists a constant C > 1 
independent of u an J n such that for any integer k, < k < n + 2, we have 

(5.4) i ( Y. < ( E ii^^"^'"'a?fiiu.)' 

|Q|=fc ^=0 |a|=£ 

^ II' '-'x "llo;>Vs J 

|a|<l 

(ii) Let m > 1 be an integer such that m + 1 > —(3. Let n > m — 1 be another integer. 
7/u G J™(We) anJ f G J^+gl^e)' ^^^^'^ " ^ JJ^'^^^Cl^') anJ ^/zere a constant C >l 
independent ofu and n such that for any integer k, m < k < n + 2, we have 

<5.5) ^(Ell'-"*'°'«:w)^£C'«{ g i(Ell'-'"*'a?<w.)^ 

|a|=A; £=m-l |a|=£ 

a|=m 

Proo/ Like in the case of Theorems 3.2 and 3.7, the proof relies on a locally finite dyadic 
covering of W and W^. The reference domains are now 

V = {xxG/C : i<|x^|<l}x(-i,l) 
V' = {x_L G /C : i - e < |x_l| < 1 + e} X (-i - e, I + e) 

and for /X G N and G Z: 

V^,. = 2-''(V+(0,0,f)) and V^,. = 2-^(V' + (0, 0, f )). 

We check immediately that 

^=1) U ^^^^ >^^^U U ^/^.- 

and that these coverings are locally finite. An a priori estimate between V^^^, and V/j,,^ is 
deduced from a reference a priori estimate between V and V' by the change of variables 
X — X = 2~^(x + (0, 0, |)) that maps V onto V/^^i, and V onto V/l,i,. Here we use the fact 
that the operators L, T and D are homogeneous with constant coefficients. Then the rest of 
the proof goes exactly as in the case of the plane sectors. □ 

5.2. Tangential regularity along the edge (homogeneous norms). The result in the previ- 
ous sections only rely on the ellipticity of the boundary value problem under consideration. 
Now we will require a stronger condition, which is a local Peetre-type a priori estimate in 
an edge neighborhood. From this condition we will derive analytic type estimates for all 
derivatives d^,^^ in the direction of the edge. 
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Assumption 5.3. Let /3 G M. Let W and W = We be the domains defined in (5.3) for some 
e > 0. We assume that the following a priori estimate holds for problem (5.2) on V = W: 
There is a constant C such that any 

u G Kim , 

solution of problem (5.2) with f G K°^2(^')' satisfies: 



(5-6) II"IIk^(w) - ^[^\^hl^,{w') + II"IIki+,(w') 

Remark 5.4. (i) Assumption 5.3 is independent of e (although the constant C depends on 
it), and more generally independent of the choice of the domains W and W, if they satisfy 
the following conditions: There exists a ball with center on the edge e contained in W, and 
W contains WnV. 

(ii) The inequality (5.6) is a Peetre-type estimate, since K^(>V) is compactly embedded in 
( in) As a consequence of Theorem 5.2, it is equivalent to postulate the estimate 

II"IIki(w) - ^(llfllK«+,{w') + II"IIki+,(w') 

for all u G Kj(>V') n Hf„,(W' \ e) 

Remark 5.5. Assumption 5.3 can be characterized by a condition on the partial Fourier 
symbol of L along the edge, as follows. If we write the system L in the form L{dy,^, dx^), 
its partial Fourier symbol L{E,) is defined on the sector K, for all ,^ G M by 

(5.7) L(0(9xJ = L(c',^,zO, x^G/C. 

We define Tj and Di on the same way. We also need the weighted spaces on /C 

(5.8) E™(/C) = {iL G ll,{K) : max{r^+l"l,r^+'"}9> G l^JC), Va, |a| < m}. 
Then Assumption 5.3 holds if (and only if) the problem 

L(0 u = f in /C 

(5.9) { T,(Ou = on F,, z = 1,2, 

A(0" = on F„ z = 1,2, 

defines a injective operator with closed range from E^(/C) into E°_^2(^) for C = il- I^i 
[34], Maz'ya and Plamenevskii introduced the spaces E^(/C) and proved that isomorphism 
properties of the transversal problem (5.9) are necessary and sufficient for Fredholm prop- 
erties of the boundary value problem (5.2) on the wedge. The same technique proves that 
left invertibility of the transversal problem (5.9) implies the semi-Fredholm estimate (5.6) 
of Assumption 5.3. 

The first step for higher order estimates is the p-estimate for which we control the depen- 
dence of the constant C in (5.6) on the "distance" between W and W. 
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Lemma 5.6. Under Assumption 5.3, let R G [0,e) and p G (0,£: — R]. Assume that u G 
K^(>V/j) zi' a solution of problem (5.2) wzY/z f G K°_,_2(V)/or V = Wi?+p. There exists a 
constant C independent of u, R and p such that 

(5.10) II"IIk2(w«) - ^(iI''IIko^,(w«+,) + ^ 'II"IIki+,(w«+p) + ^ '""IIko+,(w«+,))- 

Proof. We introduce a special family of cut-off functions Xp- Let x G ^°°(]R) be such that 
X = 1 on (— oo, 0] and x = on [1, +oo). Define Xp on M by: 

(5.11) XpW = xf^ 

Thus Xp equals 1 in [— 1 — i?, 1 + i?] and outside ( — 1 — R — p,l + R + p). Then we set 

(5.12) Xp(x) = Xp(|x±|)Xp(a;3)- 
Thus by construction, cf. (5.3) 

Xp = l on Wr and Xp = outside W/j+p. 
We note the following important bound on the derivatives of Xp 

(5.13) 3D > 0, Vp > 0, Va, |a| < 2, \d^Xp\ < Dp'^^'K 

Then in order to prove (5.10), it suffices to apply estimate (5.6) to XpU and to check that the 
commutator [L, Xp] applied to u satisfies 

(5.14) ||[L,Xp]u|| ^ < C'(p~1u|Iki (w ) + P^^\Mko (w ^|■ 
The latter estimate is an obvious consequence of (5.13) and the fact that 

for all a, |a| < 1. □ 

Corollary 5.7. Under the assumptions of Lemma 5.6, ifdxj G K°^2(^if+p)' ^^^^'^ ^^gU ^ 
K^()A^/j) and there exists a constant C > 1 independent of R, p and u such that 



Proof. For any h < p/2, we apply (5.10) in 'Wr+p/2 to mh defined by 

: X -)> /i~^(u(x + /ics) - u(x)), 
where 63 = (0, 0, 1). This yields 

(5.16) 

— 1 II II —2 II II 

where C is the positive constant from Lemma 5.6. By noticing that 



Vh = h W dx^u{x + tea) dt, 
Jo 
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we check that for all /i < p/2 

This shows that the right-hand side of (5.16) is bounded uniformly in h. Therefore passing 
to the limit in (5.16), we find that d^^u belongs to K^(>Vij) and that (5.15) holds. □ 

Corollary 5.8. Under Assumption 5.3, let u G K^(>Ve) be a solution of (5.2). Let R E 

[0,e/2] and R' > e/2 with R + R' < e. Then there exists a constant C independent of R, 
R! and u such that for all £ G N, we have 

i=o 

Proof. If £ = 0, this is a consequence of estimate (5.6). For £ > 1 the proof is divided into 
two steps. To keep notations simpler we take i? = 0. 

(i) We first prove by induction on i that if p < s/{2t — 1), then 



+ p-'\W\ 



-l)p)^^ """k^ + i(W(2£-1)p) }' 



where C > 1 is the constant from Corollary 5.7. 

• If £ = 1, the estimate (5.18) is nothing else than (5.15). Hence it suffices to show that if 
(5.18) holds for it holds for ^ + 1. 

• For that purpose, we first apply (5.18) to defined as before by 

: X /i~^(u(x + /les) - u(x)), 
and passing to the limit in h, we get 



(>V(2£+l_j)p) 



For the second term of this right-hand side, we apply (5.15) to u but between W2ip and 
W'r2£+i~)o, while for the third term we use the fact that ||9a,,u||.,i ,,,, ^ < ||u||,,2m,, ^. This 
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leads to 

IIC'"IIk|,w, s (2C)'i:''-"-^'ll«"i"llKj„,„,„,..,„, 

By the change of index j' = j + 1 in the sum on j, we finally get (since {2CY < 2^C^~^^) 

i+i 

Since C >1,C +1 <2C, and this proves that (5.18) holds for i + 
( ii) Now we choose p such that 

>V(2^_i)p C We' with £ = e/2. 

This holds if we take 

p = - With 7 = min|-,l|. 
Hence applying (5.18) with this choice of p, we obtain for alH > 1 



(5.19) " j=i 

+ (7-'mu|l,j,„.„+(7-)-'^-'l|u||,.^_,„,„}. 

Since 7 < 1, (7~^)(^~-') < (7^^)^. Moreover by Stirling's formula, one has 

/ < ! 

for some > 1. We find 



since j\ < < P. Inserting this into (5.19) gives, with Ci = 2C'y~^S, 



Using the trivial inequality £ < 2^, we arrive at 



which, combined with (5.10) between Wg/ and We , yields the requested estimate. □ 
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5.3. Anisotropic estimates in dihedral domains (homogeneous norms). We are now 

ready to prove the main results of this section, namely the weighted anisotropic regularity 
of solutions of our local boundary value problem (5.2). For this we introduce the following 
new class of weighted spaces: 

Definition 5.9. Let (3 he a real number and let m > be an integer. 

Let W be a subdomain of the dihedral domain V. We recall that r = \x± \ denotes 
the distance to the edge e = {x^ = 0}. The anisotropic weighted space with 
homogeneous norm M™(>V) is defined by 

(5.20) M'^iW) = {ue Lf„,(W) : r^+l"^l9> G l^iW), Va, |a| < m] 

where for a = (ai, ol^-, as), oli_ = (ai, 0:2) is the component of a in the direction 
perpendicular to the edge e. The norm of this space is defined as 

m 

(5-21) lkllLn.n = y y Ik^'^'^^'^x «lln • 

V / II llM™(>V) / J / X ll0;W 

fc=0|a|=fc 

Theorem 5.10. Let /3 G M. Under Assumption 5.3, let u G K^(We) he a solution of 
problem (5.2). Tff G M^^gl^e)' ^^^^ ^ ^ M^(W), and there exists a positive constant C 
independent of u and n such that for all integer k, < k < n we have 

(5.22) i(Eii'''^""^'a?"ii:„)*<c-'{t^(Eii'-'"^""^'s?fii:vv. ' 

\a\=k i=0 \a\=£ 

+ II"IIk1+,(W,) 

Proof, (i) We first apply the isotropic estimate (5.4) between W and >Ve/4, and combine 
with (5.6) between We/4 and Ws/2 (cf Remark 5.4 (i)). This yields the estimate for all k, 
< k < n, and with e' = e/2 

^ I X ^ II aj^\n\ II 2 \2 , „j.ii fx ^ J- / V ^ ., n I o I I „,| „„, _ II 2 



( E ii^^^'-'^^ii^w)' <^^^HEa5: 



\a\=k i=0 ' \a\=e 



hand sic 

/c = or 1), we write the previous inequality in the slightly weaker form 

1. ( ^ ii,.*+H<¥uii^^„)' < cr{t^ ( E ii-^*^*'a?fii:„,., 

\a\=k e=0 \a\=e 



In order to absorb the term Hr^+^f || ^ in the sum on the right-hand side (including when 
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We reduce the left-hand side to any a = {a_i_, 0) of length g > 0, and bound r^+^+l°l by 
^/3+2+|q±| jj^ ^j^g right-hand side (recall that r is bounded in Wg) to obtain for all g, < g < n 



l"±l=? 

(5.23) 



=0 \a\=£ 



(ii) We now prove that for all yU = 0, . . . , n and for all g = 0, . . . , n — /x one has the following 
estimates with k := q + ji and a constant C3 independent of u, g and jj, 



lax|=ij 

(5-24) , 

£=0 ■ \a\=i 

1. If yU = 0, this estimate is a consequence of (5.23) since We' C Wg. 

2. If yu > (or equivalently g < A;), we apply (5.23) to d'^^u to obtain 



^ I 

(5.25) 



=0 \a\-- 

The last term of this right-hand side is now estimated with the help of Corollary 5.8. Using 
that 

II^^3"IIk1+i(W,,) - ^"IIk2(W,,)' 

and applying (5.17) between We' and W^ with £ = /i — 1, we obtain 

II''-"IIk;,«(w„) s^S't*'- 1)' (i:^lia:3fllKS„,w., + MIk;,,w.,)- 

i=o 

Using this estimate in (5.25) we obtain that 

1 ( ^ ||/+i-i9-aj3u|i;„)' < cf x: I ( E ll'-''"^'-'<a?fli:.w. )' 



|Qx|=g ^=0 \a\=e. 

Jl-1 



+crcr(/.-i)!(Ej!ii^^3fiiK0^,(w.) + ii"iiK^_ 



.1II^^-3'IIk0^2(W,) ^ ll"llKi + ,(W,) 

i=o 
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Multiplying this estimate by g!(A;!) ^, we find (since — ^ < 1) 



2 



/ 1 

3=0 



Vox the first term of this right-hand side we finally notice that d^.p'^ = and that 

la + (0, 0, /i) I = £ + /i. Hence we have to check that 

q\ ^ 1 



m - (£ + /i)!' 

which is equivalent to 

m - ' 

and holds since i + fi < k and q < k. 

Altogether we have proved that (5.24) holds for all G N such that q + fi = k. 

(Hi) Summing the square of this estimate (5.24) on g = 0, . . . , /c and /i = 0, . . . , /c — g, we 
arrive at 



\a\=k 1=0 ' \a[ 



+i|"|Ik^,.(w.) 



This proves the theorem. □ 

5.4. Anisotropic estimates in dihedral domains (non-homogeneous norms). In this last 
part of section 5 devoted to local estimates in dihedral domains, we investigate the situation 
where the a priori estimate holds in the J-weighted scale instead the K scale. We set: 

Assumption 5.11. Let /3 G M. Let m > 1 be an integer such that m + 1 > —fi. We assume 
that the following a priori estimate holds for problem (5.2): There is a constant C such that 
any 

UG J^+HW), 

solution of problem (5.2) in V = W with f G J!^!r2^(VV"), satisfies: 



(5.26) ll"llj™+Vw) ll^llj--(w') + ll" 



'(w) - ^ VII- ||JJ+V(W') ' ii"iij^+i(w'). 
Remark 5.12. Using the analogue of Proposition 3.5 for dihedral domains, we obtain that in 

im 



the situation of Assumption 5. 11 the norm in the space S^'^^ ( W) is equivalent to 



(5.27) ( ||r'"'"^''+l"l'°>9>"' 



0; W 

|a|<m,+l 
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Remark 5.13. Along the same lines as Remark 5.5, we have a characterization of Assump- 
tion 5.1 1 by the partial Fourier symbol of L: Assumption 5.11 holds if (and only if) problem 
(5.9) defines a injective operator with closed range from J™"*"^ (/C) into JJ^2 (^) for ^ = ±1 . 

The non-homogeneous anisotropic weighted spaces are defined as follows on the model 
of the homogeneous ones (Definition 5.9): 

Definition 5.14. Let /3 G M. Let n > 1 be a natural number such that n > — /5. 

Let W be a subdomain of the dihedral domain V. The anisotropic weighted space 
with non-homogeneous norm N^(VV) is defined by 

(5.28) N;^(W) = {ue : r'^--{^~^\»^lo}gy ^ L^{W), Va, |a| < n} 

endowed with its natural norm. 

Our aim is to prove the "non-homogeneous" analogue of Theorem 5.10: 

Theorem 5.15. Let /3 G M. Let m > 1 be an integer such that m + 1 > — /3. Under 
Assumption 5.11, let u G J^iWe) be a solution of problem (5.2). T/'f G H'p^2(^e) for an 
integer n > m, then u G N^(>V), and there exists a positive constant C independent ofu 
and n such that for all integer k, < k < n we have 

\a\=k 

k 1 

E- ( \^ II max{/3+2+|ax|,0}oar||2 \2 || ^ 

Proof We first notice that Theorem 5.2 yields u G ^'^'^{We') for any e' G (0,£:). As 
n > m — 1, we have obtained the basic regularity u G JJ~^^{yVe')- We review now the 
sequence of steps leading to Theorem 5.10 and adapt them to non-homogeneous norms. 
(i) Applying (5.26) to XpU with the function Xp introduced in (5.12), we obtain, - compare 
with (5.10), 



lull .m+l. 



m 
A=0 

(ii) By the differential quotients technique we deduce, - compare with (5.15), 

m 

A — 

since H^a^guH ^ ^ is bounded by {|u{| ,,„+i-a,,., v 



'/3 + l+A 
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(Hi) Iterating this on the model of (5.18) we find for > 1 

m 
A=0 

leading to the analytic type estimate, - compare with (5.17), 

(5.30) 77l|9f3"ll r"+i(w ) ^^'^'{E^ll^-3^"llj'"-irw ^ + \\"h^ rw )]■ 



f/vj To prove (5.29), we start with the proof of, - compare with (5.23) 

^ \ I 

ax =9 

(5-31) , 

E— f \^ II max{/3+2+|Qxl,0}nar||2 ^ ,, ,, 1 

€=0 |Q|=f 



• For g = 0, . . . , m, we rely on the estimate (5.26) combined with the use of the norm (5.27) 
for J^~'~^(>V): If we restrict the left-hand side to the derivatives of the form d^-^ and replace 
the weight r"''^^{l^+^+H:0} ^max{/3+2+|axl,o} in the right-hand side, we obtain (5.31). 

• For q > m + 1, we combine the estimate (5.26) with the isotropic non-homogeneous 
estimate (5.5) and making the same restriction to d^^ in the left-hand side and the same 
change of weights in the right-hand side. 

(v) We continue with the proof that for all /i = 0, . . . , n and for all g = 0, . . . , n — one has 
the following estimates with k := q + n and a constant C3 independent of u, g and jj, 



\a±\=q 

(5.32) , 



— { \^ ||^max{/3+|axl,0}QaxAM 1,11^ V < (7 



E — { \^ II max{/3+2+|ax|,0}aar||2 

e=o ■ \a\=e 

1. If /i = 0, this estimate is a consequence of (5.31) since We' C We. 

2. If yU > (or equivalently q < k), we apply (5.31) to df^^u to obtain 

1 / ^ 9 • ' 



E— f \^ II max{,3+2+|ai|,0}oao^ r||2 \ 2 , 



\a± =q 

(5.33) , ^ 



=0 |al=£ 
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The last term of this right-hand side is now estimated with the help of (5.30) with i = — l 
Using this estimate in (5.33) we obtain that 

— ( ||^max{/3+|axl,0} aax/5M 1,11^ V< 

\o:±\=g 

e=o ' \a\=e 

+ Crc^i^ - 1)! {Y.lKJhj;H„., + )■ 

j=0 

We note that the norm in the space J^^2 0^£) equivalent to (cf. (5.27)) 

Q|<m.— 1 

Thus dividing the latter estimate by /i! and recalling that k = q + fx we deduce 

\a±\=q 

£=0 \a\=e 

_Lr;k+lf\^±_ II max{,3+2+|a|,0}aQQi fll _l||„|| ^ 

+ ^5 I^Z^ -1 ll'^ ^=<^^3^ll0;W, + ll"llJ^+l(W,)y■ 

j=0 ■ |a|<m-l 

From this we deduce (5.32). The final way to (5.29) is very similar to the conclusion of the 
proof of Theorem 5.10. This ends the proof of Theorem 5.15. □ 

Remark 5.16. We note some similarities between our estimates and those obtained in [24] 
for the Laplace operator. Our argument based on the dyadic partition technique clearly 
improves the structure of the whole proof. 

6. Natural anisotropic weighted regularity shift in polyhedra 

6.1. Edge and corner neighborhoods. Let f2 be a polyhedron in R^, that is a domain 
whose boundary is a finite union of plane domains (the faces Fg, s G y). The faces are 
polygonal, the segments forming their boundaries are the edges e of fi, and the ends of the 
edges are the corners c of ^2. We denote the set of edges by and the set of corners by ^. 
Edge openings may be equal to 27r, allowing domains with crack surfaces. 

In order to prove global regularity results in suitable weighted Sobolev spaces, we intro- 
duce comer, edge and edge-vertex neighborhoods of Cl. For a fixed comer c G ^, we denote 
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by S'c the set of edges that have c as extremities. Similarly for a fixed edge e G c?, we denote 
by the set of corners that are extremities of e. Now we introduce the following distances: 



(6.1) 



rc(x) = dist(x, c) 



(x) = dist(x,e), pce(x) 



Tp X 



rrX 



There exists e > small enough such that if we set 



(6.2a) 



{xen 
{xen 



re(x) < E and r^x) > e/2 Vc G K}, 
rc(x) < e and pce(x) > e/2 Ve G 4}, 
rc(x) < e and pce(x) < e}, 



we have the following properties: 



(6.2b) 



B{c,e)nB{c',e) 



0, Ve' 7^ e, 
0, Vc' ^ c, 
0, VeVe. 



We also define the larger neighborhoods with e" < e < e' 



(6.2c) 



{xen 
{xen 
{xen 



re(x) < e' and r^x) > e" /2 Vc G %}, 
rc(x) < e' and pce(x) > e"/2 Ve G 4}, 
rc(x) < e' and Pce(x) <£:'}, 



assuming the e' and e" are sufficiently close to e for the above properties (6.2b) to hold 
for fig, n'^, and fi'^g. We finally introduce the smaller neighborhoods fig, fi^', and fi^g by 
inverting the roles of e' and e" and set. 



(6.2d) 



fi 



fi. 



fi^. 



U U 



We define fin as the remainder: 



(6.2e) 



fio = fi \ fi<:^ u fi^ u n^s 



Note that fio is far from the singular points of fi. We finally choose a larger "smooth" 
neighborhood fi[) C fi such that fi n fio C fig and fi[) n ((^' U ^) = 0. 

Let V be any subdomain of fi. We consider the system of local interior and boundary 
equations 



(6.3) 



f in fin V, 

on n V, s G ^, 
on n V, s G ^, 



where the operators L, Tg and are homogeneous with constant coefficients and form an 
elliptic system. The choice V = fi gives back the global boundary value problem on the 
polyhedron fi. 
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Definition 6.1. OnV CQ, for m eN and § = {/^cjce-ir U {Pejess-, the weighted space with 
homogeneous norm K™(V) is defined as follows, cf. [37, 38, 12, 13] 

(6.4) K|^(V) = |m e Lf„,(V) : Va, |a| < m, 9> G L2(V n fio) and 

rc(x)'^^+i"i a> el\vn Qc) Vc e ^, 
re(x)'^^+l"l 9> G L2(V n Qe) Ve e , 

rc(x)'''+l"l Pce(x)^^+I"l dy G L2(V n l^ce) Vc G ^, Ve G 4}, 
and endowed with its natural semi-norms and norm. 

Note that the condition in the edge-vertex neighborhood Vice can be equivalently written 



as 



rr X 



,/9c-/3e 



Remark 6.2. The semi-norms issued from (6.4) are equivalent to the globally defined semi- 
norms 

{Ell{n'-c**'}{n(^)''""'V>tj*. k=o,....n,. 

\a\=k ceV eeS 

Here r^^ denotes the distance function to the set of corners. With this expression, the 
relations between our spaces K^(f2) and the spaces V^f{^) defined in [37, §1.2] or [38, 
§7.3] become obvious: 

(6.6) K^in) = V^fin) if p = 2, /3={/3, + m}^^.^, 5 = {/3e + m}^^^ . 

6.2. Anisotropic weighted spaces with homogeneous norms. Unlike in the conical case, 
the weighted spaces K™ are in a certain sense too large to describe accurately the regular- 
ity of solutions of the elliptic problem (6.3) along the directions of edges. Mimicking the 
definition of the spaces M™ in the pure edge case, cf. (5.20), we particularize for each edge 
e E S', the derivatives in the directions transverse or parallel to that edge by the notations 

(6.7) dy (transverse) and a^n (parallel), (e G f?), 
so that 

= dy . 

Of course these directions are edge dependent. They are well-defined in each of the domains 
and Vice determined by the edge e. 

The following spaces were introduced in [12, 13] for similar purposes: 
Definition 6.3. On V C 1^, for m G N and j3 = {Pcjce'^f U {(3e}ee<^, we define 

(6.8) M™(V) = |m G ll^iV) : Va, |a| < m, dy G L\V n Qq) and 

rc(x)^^+i"i dy G i\v n Qc) Vc g 

re(x)^^+l"^l dy G l\V n Qe) Ve G 

rc(x)^^+l"l Pce(x)^^+I"^l dy G l\V n fice) Vc G ^, Ve G C^c], 
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We denote by 11 ■ |L, „ and I ■ L , its norm and semi-norm, namely 

^ II llM;m,^;V I lM;m,^;V ' ^ 

m 



iM;m,aV I iM;£,aV 

i=0 



with 



(6-9) l-lL;,,;V = E (ll^>llo;Vn.o + E ll-^"'"'^x «llo;Vn.. 

\a\=l cG<r 

Note that the condition in the edge- vertex neighborhood fice can be written equivalently as 
We can then define the corresponding analytic class as follows: 



Definition 6.4. We say that u G A^(r2) if u G Mo(r2) for all /c > and there exists a positive 



'II 

constant C such that 

I lM;fc,^;n — — 

We rephrase Assumption 5.3 for the dihedral neighborhood Vt^: 

Assumption 6.5. Let e G and (3e G M. We assume the following a priori estimate: There 
is a constant C such that any 

UGKj^(fie), 

solution of problem (6.3) in V = with f G K°^_^2(^e)' satisfies: 



We can apply Theorem 5.10 to the edge neighborhood Vt^- We obtain that under As- 
sumption 6.5, any solution u G K^^(r2g) of problem (6.3) with f G ^'p,+2i^'e) satisfies the 
uniform estimates for < A; < n 

\a\=k i=<d \a\=l 

Now we consider the edge- vertex domain Vt^e- 

Proposition 6.6. Let c G ^ and e G S^. Let 13 = {/3c, /3e}- Under Assumption 6.5, 
any solution u G K^(r2ce) of problem (6.3) wzY/z f G M^^_2(fice) belongs to M^(f2ce) and 
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Figure 2. Nested edge neighborhoods (section determined by azimuthal 
angle 9^ = constant) 



satisfies the uniform estimates for < k < n 

(6-12) ^ ( E ||r.^-l-lp^-l-l9>||^^,J^ < C''-^{ 

\a\=k 

E ^ ( E ii'-^^*'°'pa«^'-'a?fii:„, )* + iiuii,.^,„, , }. 

£=0 ■ \a\=i 

Proof. We mimic the proof of Theorem 3.2. The proof of estimate (6.12) is based on a 
locally finite dyadic covering of Qce and D.'^^. Define, compare with (6.2a)-(6.2c), 

V = {x G : I < rc(x) < £ and pce < 
V' = {xen : ^ < rc(x) < e' and pce < ^'j, 



and for G N: 



We check: 



Vf, = 2-^'V and Vl = '2'^V. 



nce = \JV^ and 1];,= U V;. 



The estimate (6.1 1) between fie and also holds in the configuration of V and V' which is 
similar: V and V' are nested edge neighborhoods which do not touch any corner, see Fig. 2. 
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Since Tc is bounded from above and from below by strictly positive constants, the distance 
Te is equivalent to pce on the reference domains: We have 

\a\=k e=0 \a\=e 

la|<l 

for any reference function u satisfying the boundary conditions of (6.3) and f := Lu. 

For the same reason, we can insert powers of Tc in the above estimate, to obtain our new 
reference estimate 

\a\=k 

£=0 ' \a\=i 

+ 5^ ||rc(x)^^+Hp^^(x)/^e+H5.^||^^„J. 

|a|<l 

The change of variables x — )• x = 2^^x maps V to (resp. V' to V/1). We note that 

Pce(x) = pce(x) and rc(x) = 2^rc(x). 
With the change of functions 

u(x) := u(x) and f(x) := Lu, which implies f(x) = 2~^^f(x), 
we deduce from estimate (6.13) that 

\a\=k 

^i. 2M(^c+2)(^ ^ 2-^'^||r.(x)^+^+l-lp„(x)^^+^+l--l9,"f||^^^,)^ 
e=o ' \a\=e 

+ (2^^^ Yl Ikc(x)^^+I"lpce(x)^^+I"l9>||^^^,)'}. 

|q|<1 

Multiplying this identity by 2^^^, taking squares, and summing up over all /i, we get the 
requested estimate (6.12). □ 

The estimates in pure vertex domains fic (i-C, close to comers but "relatively far" from 
the edges) are similar to those in obtained in Theorem 3.2 for plane sectors: 
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Proposition 6.7. Let c G ^ and § = {/3c}- Any solution u G K^(^^c) of problem (6.3) with 
f G M^;^2(^c) belongs to M^(r2c) and satisfies the uniform estimates for < k < n 

\a\=k 

Proof. The proof is again based on the argument of dyadic partitions with reference domains 
defined as 

V = {xGfic, f<rc(x)<£} and V' = {x G fi^, ^<r,{x)<e'}. 
and for G N: 

V;, = 2"^V and = 

We check: 

= U ai^d = U 

We can apply the a priori estimates of the smooth case between V and V, cf. (3.9) and 
deduce (6.14) in the same way. □ 

We obtain now the anisotropic regularity shift in homogeneous weighted spaces on poly- 
hedra : 

Theorem 6.8. Let Q be a polyhedron and (3 = {/3c, Pe} be a weight multi-exponent. Let 
Assumption 6.5 be satisfied for all edges e E Let u G Hfo^(f2 \ S') be a solution of 
problem (6.3) in V = Vt. Then the following implications hold 

(6.15a) u G K\{n) and f G M™_^2(^) =^ " ^ M^in) (m G N), 

(6.15b) u G K^(fi) and f G A^+2(fi) ^ u G A^(f]). 

Proof. The proof is a consequence of 

(i) elliptic estimates in the smooth case applied between VLq and Vl'q, 

(ii) pure corner estimates (6.14), 

(Hi) edge estimates (6.1 1) between the pure edge domains He and fig, 
(iv) edge- vertex estimates (6.12). 

□ 

6.3. Anisotropic weighted spaces with non-homogeneous norms. For the same reason 
as in the two-dimensional case, it is valuable to have alternative statements to (6.15a) and 
(6.15b) in which the a priori condition u G K^(fi) can be replaced by the weaker condition 
u G Jj(fi). 
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Definition 6.9. For ^ = {/3c, /3e} and n G N, let us introduce the isotropic weighted space 

(6.16) J|(V) = |n G ll^{V) : Va, |a| < n, 9> G L2(V n l^o) and 

rc(x)^'+" 9> G L2(V n fic) Vc G ^, 
re(x)'^^+" 9> G L2(V n fie) Ve G ^, 

rc(x)^'+" Pce(x)^^+" 9> G L2(V n fice) Vc G ^, Ve G ^^'c } , 

and its anisotropic companion, for n > — min{mincg<r /3c, minee^ /3e}, c/ (5.28) 

(6.17) N|(V) = |m G Lf„,(V) : Va, |a| < n, <9> G L2(V n fio) and 

^^^^)max{/3c+|a|,0} ^ |_2(y ^ y^, ^ 



^^)max{/3e+|axl,0} ^ l2(V H fice) Vc G ^, Vg G ^c}- 



We note that, like in the case of K-weighted spaces, the semi-norms issued from (6.16) 
are equivalent to the globally defined semi-norms, compare with (6.5) 

2 

r i. - n 



(6.18) {j:\\{Urs-*"}{u&ni^-4j". *=0 

\a\=k eGff ^ 

It is useful to introduce, in the same spirit as in [38], a full range of intermediate spaces 
between K^(fi) and J^(fi). 

Definition 6.10. Let us flag a subset of comers and a subset (oq of edges, and define 
J^(V; ^0, ^o) as the space of functions such that all semi-norms 

(6.19) ||{ n n --i'-""}! n (^;)""'°'}{ n O'""} 



0;V 



are finite for |«| < n. Anisotropic spaces N^(V;^o,A) are defined similarly, replacing 
in (6.17) the weight ^^^^{z^'+l^l'Oi by r^+'°' when c G ^o, and {re, pce}"'^''^^^+'"''°^ by 
{^6,^06}^^'+'"' when e G 4. The sum of the squares of these contributions for | a | = n 
defines the squared semi-norm 

I |2 

I In|(V;%,^o)' 

Note that with = = we obtain the maximal spaces already introduced in (6.16) 
and (6.17): 

(6.20) J|(V) = J|(V; 0, 0) ; N^(V) = N|(V; 0, 0) . 

The corresponding analytic class is defined as usual: 

Definition 6.11. We say that u G B^(fi; ^o, 4) '^i u ^ N^(fi; <^o, 4) for all k > := 
— min{mincg<r /3c, minee^f /3e} and there exists a positive constant C such that 
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In accordance with (6.20), we write B^(r2) for B^(r2; 0, 0). 

Remark 6.12. (i) Choosing = ^ and S'o = S, we find that the spaces J|(fi;^, ^), 
N^(fi; (g) and B^(l]; S) coincide with the homogeneous spaces K|(fi), M^(fi) and 
A^(fi), respectively. 

{ii) The following relations hold between our spaces J^(fi; ^o, <^o) and the spaces W^i'(fi) 
of Maz'ya and Rossmann [38]: 

(6.21) }^l{^-^,Qi) = WJI{^) ifp = 2, /3 = {/3c + m}^^,^, 5 = {/^e + m}^^^ . 

In these spaces, the non-homogeneity is only related to edges. Under the same condition 
as in (6.21), the intermediate spaces J) of [38, §7.3] coincide with our spaces 

J™(fi; (oo) if is chosen as the same set of edges as J. 

{in) Our analytic class B^(f2) coincides with the so-called countably normed spaces B^^(^l) 
introduced by Guo in [21]: If Guo's edge and comer exponents /3ij G (0, 1) and /3m G (0, ^) 
satisfy /3ij = /3e + ^ and /3m = /^c + respectively, then -B^(fi) = B^(fi). 

We state the assumption for J -weighted spaces corresponding to Assumption 5. 1 1 for the 
dihedral neighborhood f2e: 

Assumption 6.13. Let e G ^. Let /3e G M. Let m > 1 be an integer such that m + 1 > —/3e. 
We assume the following a priori estimate: There is a constant C such that any 

uG jr'(^e), 

solution of problem (6.3) in V = with f G J)^+2(^e)' satisfies: 



We then have the following anisotropic regularity shift result in the non-homogeneous 
weighted spaces N|(fi; 'rf, 0) and B^(fi; ^, 0): 

Theorem 6.14. Let Q be a polyhedron and § = {/3c, /3e} a weight multi-exponent. Let 
m > 1 be an integer such that m + 1 > —fSefor all edges. Let Assumption 6.13 be satisfied 
for all e G <f . Let u G H^Q^(r2 \ S) be a solution of problem (6.3) in V = Vt. Then the 
following implications hold 

u G JJ(r];^,0) anJ f G N|+2(^;'^.0) =^ ugN|(1];^,0) (n > m), 
u G J|'(fi; ^, 0) anJ f G B^+2(^^; ^,0) ^ u G B^(fi; '^f , 0). 

Proof. The proof is a consequence of suitable a priori estimates with analytic control in the 
four types of regions in the polyhedron: 

(i) Elliptic estimates in the smooth case can be applied between VLq and 0!^. 

(ii) Pure corner estimates (6.14) are valid here: We note that in the pure comer region VLc 
the norms in K and J spaces, or in M and N spaces, are the same. 

(Hi) The edge estimates (5.29) are valid between the pure edge domains fie and Vt'^. 
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(iv) Finally, edge- vertex estimates are proved by the dyadic partition argument starting from 
the same reference domains V and V' as in the proof of Proposition 6.6. The reference 
estimate can be written as 



_|_ ^ ||^max{/3e+|a|,0}^Q!|J| 



1 

2 \ 2 



0; V 



\a\<m 

Since rc and (rc)~^ are bounded on the reference domains, we can 

• replace by pce 

• insert powers of 

in the previous estimate, thus obtaining 



\a\=k 



£=0 ■ \a\=e 



J2 lkf^^'"'Pr^''^^'"''°^5xB|lo;V' 

|a| <m 

Owing to the homogeneity of the weights with respect to Vc, the dyadic partition argument 
yields the desired edge-vertex estimate, which allows to conclude the proof of the theorem. 

□ 

Remark 6.15. (i) If we replace Assumption 6.13 by Assumption 6.5 for edges e in the 
flagged subset do, we can prove, instead of (6.23), the implications 

u e J^(fi; ^, 4) and f G N|+2(^; ^, 4) =^ u G N^(fi; ^, 4), 
u G (fi; ^, 4) and f G B^+2(f^; ^, 4) =^ u G B^(l]; 4)- 

(ii) Under Assumption 6.13, the implications in the maximal non-homogeneous spaces, i.e., 
with ^0 = 4 = 0, are also true: 

u G J^{Q) and f G N^+^l^) =^ u G N|(fi), 
^ ■ u G J^(fi) and f G B^_+2{n) =^ u G B^(l]). 

If /9c > — I for any comer c, the statements (6.23) and (6.26) coincide, since in this case the 
spaces 0) and J™(1^) are the same (consequence of Hardy's inequality). In the 

general case (6.26) can be proved by two different methods: 

• Deduced from (6.23) by an argument of corner asymptotics (at each comer, the 
asymptotics modulo J™(fi; 0) contains only polynomials): For instance when 
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m = 1, if /3c G (— |, — |) for all comers c, any element of u G J™(^^) splits as 
u = Uc + Wc in l^c, with Uc G J|'(fi; ^; 0), Wc G C^, 

and we can apply (6.23) locally near each corner, to each function Uc- 
• Directly proved by the same method as for Theorem 6. 14, starting with the reference 
estimate for A; > m 

\a\=k 

k , 1 
i=m-l ^' \a\=e ^ 

\a\=m 

instead of (6.24): The norm present in (6.24) is replaced here by the correspond- 
ing semi-norm, cf. Corollary 2.2. 



7. Analytic weighted regularity for solutions of coercive problems 

In this section, we show how Theorems 4.4 and 4.5 in the polygonal case, or Theorems 
6.8 and 6.14 in the polyhedral case, apply to solutions of variational problems. For second 
order boundary boundary value problems that allow a coercive variational formulation, one 
knows basic regularity in weighted Sobolev spaces in a form that fits the hypotheses of 
our natural regularity shift results. For polygons, this is obtained by means of Kondrat'ev's 
classical theory, and for polyhedra, such results were proved by Maz'ya and Rossmann [38]. 
As a consequence, we obtain analytic regularity for solutions of variational problems. 

Let be a polygon or a polyhedron. In coherence with the previous sections, we consider 
a sesquilinear form a, homogeneous of order 1 and with constant coefficients acting on 
vector- valued functions with components 

N N 

(7.1) = E E E E / ^y^i"") ^x^^ w d'^' 

i=l j=l \a\=l |7| = 1 

and a subspace V of H^(f2)^ =: H^(r2) defined by essential boundary conditions on the 
sides Fs of 

(7.2) V = {u G : D,u = on F^, s G ^}. 
We assume that the form a is coercive on V: 

3c, C>0, VugV, Rea(u,u) > c||u||'.^ -C||u||^.^. 
Standard examples of such sesquilinear forms are the gradient form for scalar functions 

av(M,'y)= / Vm(x) ■ Vu(x) dx 
Jo. 
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and the stress-strain sesquilinear forms in linear elasticity: 

aeia = / o-(u)(x) : £:(v)(x)dx, 

Jo. 

where e is the symmetrized gradient tensor and a = Ae, where A is a material tensor with 
the usual symmetry and positivity properties. Variational spaces V on which ay is coercive 
can be defined by any subset of the set of sides y-. 

V = {n G : u\^^ = Vs G Yd}- 

As for aeia we can take for V any space of the type 

(7.3) V = {u e H^(l]) : u|^^ = Vs G u ■ n|^^ = Vs G 

and u X n|j^ = Vs G ^n}, 

where n is the outward unit normal vector to Fs, and =5^z), and yN are disjoint subsets 
of y. As a consequence of Kom's inequality, aeia is coercive on such spaces V. 
We consider the variational problem 

(7.4) Find ugV such that Vv G V, a(u,v)= / fvdx. 

Jn 

With L, Ts and Dg defined in an obvious way, solutions of (7.4) satisfy 

{L u = f in 
Ts u = on Ts, s G y, 
DsU = on Ts sey. 

Having the natural analytic regularity shift results of Theorems 4.4, 4.5, 6.8 and 6.14 at 
hand, the issue is to find suitable exponents [3 so that 

(1) A^(r2) or B^(f2) are compactly embedded in ((]), — in order to be useful in error 
analysis for example. 

(2) Variational solutions u with sufficiently smooth right hand sides belong to K^(f2) or 
4(fi). 

Condition (1) of compact embedding is satisfied on two- and three-dimensional domains 
for all ^ < — 1 (this means that all components (3^ and (3^ are < —1). 

Condition (2) of initial regularity is the main question discussed in the rest of this section. 

7.1. Regularity of variational solutions in polygons. Let f2 be a polygon with vertices 
c G ^. We recall that denotes a neighborhood of c satisfying (4.2a)-(4.2b), ICc is the 
infinite sector which coincides with Qc near c, and (tc, 9c) are polar coordinates centered at 
c. Finally let Gc denote the set of corresponding angles 9c ■ Denoting by yc the set of face 
indices such that c belongs to the closure of Ts, the localized version of problem (7.5) near 
the comer c is 

{L u = f in /Cc 
Ts u = on Ts, s G yc, 
DsU = on Ts s G yc- 
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The standard Sobolev space H^(fi) coincides with Jl_^(f2), see (4.4). From Remark 3.4, 
we know that for the comparison of Jii(l^) with KLi(fi) we are in a critical case, namely a 
function u E H^(r2) neither has point values at corners nor satisfies r~^u G L^(r2) in general 
(see [30]). There holds 

K\{Q) c C Kii+,(fi), We > 0. 

Taking the essential boundary conditions into account that define the variational space 

V C H^(f2), one will sometimes find that V is embedded in K\(Q). This happens in 
particular if each corner lies on at least one side on which Dirichlet conditions are imposed. 
In the general case, one will just have V C K\_^_^{n) for all e > 0. Necessary and sufficient 
conditions for the embedding V C K\(n) are discussed in [16, Ch. 14]. 

a. Case V C K\{n) (homogeneous norms). In this case the analytic regularity shift 
(4.13b) in classes A^(f2) can be applied to variational solutions with well chosen weight 
exponents (3c < —1 as we explain now. For each comer c, the optimal condition on (3c is 
related to the spectrum cr(2lc) of the "Mellin symbol" 2lc of the system (L, Tg, Ds)^ at c (see 
[29, 30]): 

(7.7) a(2tc) = {A G C, 3ip e H^(Gc), <P^0, such that 

u := r^ip{6c) solves problem (7.6) with f = on /Cc}- 
Then we define 6c a, V) as the supremum of the numbers 6 > such that 

(7.8) {A G C : < Re A < 6} n a(2lc) = 0. 

As a consequence of the coercivity of the form a on V, the number 6c(^? a, V) is positive. 

Theorem 7.1. Let Q be a polygon. We assume that the form a is coercive on V, and that 

V C K\^{yL). If the following condition holds for the exponents (3 c 

(7.9) < -/3c - 1< bc{Vl, a, V) Vc G 

then any solution u G V of the variational problem (7.4) satisfies the regularity result: 

(7.10) fe^p_+2{^) =^ ueA^in). 

Proof. Invoking the general theory of comer problems in the variational setting, we know 
that if (7.9) holds, then 

(7.11) feK°+2(fi) =^ uGK|(fi). 

The proof of this essentially goes back to Kondrat'ev [29], see also [16, Ch. 10] for more 
details on the application of Kondrat'ev's technique to variational problems. Then (7.10) is 
a consequence of (4 . 1 3b) and (7.11). □ 

Example 7.2. Let us consider the gradient form a = ay on scalar functions. The associated 
operator is the Laplacian A. Let ooc be the opening of ^2 near the vertex c and denote by F* , 
2 = 1,2, the two sides of Vt containing c. 



'2lc is also called "operator pencil" generated by the system (L, Tg, D^). 
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(i) For the Dirichlet problem, we have V C Ki^(f2) and 

6e(l^,av,H^) = -, CG^. 

UJc 

( ii) In the mixed Neumann-Dirichlet case, if at all corners Dirichlet conditions are imposed 
on at least one side containing c, we still have V C K\{n) and 

6c(f^,av, V) = — if c G and 6c(^, av, V) = if c G ^m, 

where is the set of Dirichlet comers c (Dirichlet conditions on both sides F*) and 
the set of "Mixed" corners c (Dirichlet conditions on only one side F*). 

b. General case (non-homogeneous norms). If V ^ K\^{Vl) or for more general data with 
a nonzero Taylor expansion at corners, it is advantageous to use the analytic regularity shift 
(4.15b) in classes B^(fi). Let us recall from formula (4.10) that for -/3c - 1 G (A;, A; + 1) 
(with a natural number k): 

Theorem 7.3. Let Q be a polygon. We assume that the form a is coercive on V. If condition 
(7.9) holds for the exponents f3c, then any solution u G V of the variational problem (7.4) 
satisfies the regularity result: 

(7.12) f G B^+2(fi) ^ u G B^(fi). 

Proof. The proof relies on regularity results in spaces with non-homogeneous norms: By 
a modification of Kondrat'ev's method, see [35, 30] and [19], one can prove that if (7.9) 
holds, for any m > max{— /3c} we have the implication 

(7.13) fGJ^-/(^^) =^ ueJ^in) 

for variational solutions: In addition to the standard theory, polynomial right-hand sides of 
degree [— (3^ — I] — 2 at each comer c have to be taken into account. In dimension two of 
space, the condition that the problem (7.6) with a polynomial f of degree [—(3c — 1] — 2 
has a polynomial solution u on the infinite cone /Cc is a consequence of the condition < 
—/3c — 1 < 6c(f2, a, V). A complete proof in this framework is presented in [16, Ch. 13 & 
14]. Then(7.12)isaconsequenceof (4.15b), and (7.13). □ 

Example 7.4. Let us come back to the gradient form a = ay on scalar functions. For any 
mixed Neumann-Dirichlet problem, including the pure Neumann problem. Theorem 7.3 is 
valid and we find 

bcin, av, V) = — if c G "Tfl U and bdCl, ay, V) = if c G ^m, 

UJc 2Cc'c 

where '^^d is the set of Dirichlet comers, is the set of Neumann corners, and the 
set of "Mixed" comers c. Thus 6c (^7 ay, V) will always be greater than |. For the pure 
Dirichlet or pure Neumann problem on a convex polygon, it will be greater than 1, and for 
some triangles even greater than 2, but never greater than 3. 



47 



Remark 7.5. Theorem 7.3 has to be compared with earlier results by Babuska and Guo'': The 
Laplace operator with non-homogeneous mixed boundary conditions is considered in [4, 5]; 
more general scalar second order operators with analytic coefficients are addressed in [3], 
and finally the Lame system of linear elasticity with non-homogeneous mixed Dirichlet- 
Neumann boundary conditions is investigated in [22]. These results are at the same time 
more general than Theorem 7.3 since they address non-homogeneous boundary conditions 
and variable coefficients, but more restrictive since they do not include a full class of coer- 
cive second order systems with a unified approach. 

7.2. Regularity of variational solutions in polyhedra. Let f2 be a polyhedron with edges 
e E and corners c G ^. We recall from (6.2a)-(6.2b) the edge neighborhoods and the 
corner neighborhoods fl B{c, e). Then /Cc is the infinite cone which coincides with Q in 
B{c, e), and Gc denotes the set of corresponding solid angles 6^ = {x — c)r^^. For any edge 
e, let We be the wedge coinciding with in fie and JCe be the plane sector such that such 
that We = /Ce X M. 

The comparison between the variational space V and weighted spaces K^(fi) and J^(fi), 
cf. (6.4) and (6.16), still involves the multi-exponent I3c = [5e = —1 and essential boundary 
conditions: We have 

and, in the Dirichlet case 

Hj(fi) c K\{n). 

Moreover, the intermediate space 

^, 0) = {ue : r~^u G Vc G 

also coincides with H^(f2) by virtue of Hardy's inequality in three-dimensional cones. 

If we want to establish that weighted analytic regularity results hold in polyhedra, we 
have two tasks: 

(1) Verify Assumptions 6.5 or 6.13, which are closed range properties along the edges, 

(2) Give conditions for variational solutions to belong to spaces K^(fi) or J^(fi). 

As a matter of fact, the condition which ensures the regularity of variational solutions 
implies Assumptions 6.5 or 6.13. Hence we focus on conditions for the regularity. There 
are not so many results on regularity for elliptic boundary value problems in polyhedra. Let 
us quote [32, 33] for early results in general n-dimensional polyhedral domains in spaces 
of K type, [19] in n-dimensional polyhedral domains in standard Sobolev spaces, and more 
recently [38] in 3-dimensional polyhedral domains in spaces J^'('^, S'q), cf. Remark 6.12 (ii). 

The latter results, especially [38, Thms.7.1 &7.2], fit exacdy our requirements, namely 
in the form (6.25). For this reason we formulate Theorem 7.8 in the somewhat restricted 
framework of [38]: that is mixed Dirichlet-Neumann boundary conditions for second order 
systems. 

The regularity conditions depend on the position of the spectra cr(2lc) and cr(2le) of the 
Mellin symbols 2lc and 2le of the system (L, Tg, Ds) at the corners c and the edges e, respec- 
tively. The set cr(2lc) is defined by (7.7) on the three-dimensional cone /Cc- The set o"(2le) 



^When -/3 G (1,2), our space B^(0) coincides with their space iJ^^jC^)' 
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is defined similarly by problem (5.9) for = posed on the plane sector /Ce, namely as the 
set of exponents A G C such that the totally homogeneous problem on fCe has a nontrivial 
solution homogeneous of degree A. 

Definition 7.6. Let cr(2lc) and (T(2le) denote the spectrum of the Mellin symbol 2tc and 2le 
of the system (L, Ts, Ds) at the corner c and the edge e, respectively. Then for any edge e 
we define he{^, a, V) as the supremum of the numbers 6 > such that 

{A G C : < Re A < 6} n ^(Sle) = 

and, for any comer c, 6c(^, a, V) as the supremum of the numbers h such that 

{A G C : -| < Re A < 6} n (7(2lc) = 0. 

Remark 1.1. It is a consequence of the coercivity of the form a that the numbers 6e(^, a, V) 
are positive and 6c(fi, a, V) > — |. 

With these notations at hand, we can state 

Theorem 7.8. We consider a mixed Dirichlet-Neumann problem (7.4), which means that 
the variational space is of the form 

V = {u G H\Vl) : u|^^ = 0, S G yo}. 

Let be the set of edges e which are the sides of faces Fg with s G S^d- We assume that the 
form a (l.\) is coercive on V. If the following condition holds for the exponents f3e and /3c 

'o < -/3e - 1< &e(fi,a,V) VeG^ 



(7.14) 



-i < -/3c- 1 < &c(l^,a,V) VcG^ 



2 — 2 

then any solution of the variational problem (7.4) satisfies the regularity result: 

(7.15) f G B^+2(l^;^,4) =^ uGB^(fi;'<f,A). 

Proof First, the Fredholm Theorem 7.2 of [38] guarantees that the Assumptions 6.5 (if e G 
S'o) and 6.13 (if e G i^' \ f?o) are satisfied for any /3e satisfying (7.14). Second, the regularity 
Theorem 7.1 of [38] shows for any m > max{— /3e} the regularity u G J™(f2; ^, <S'q) with 
§ satisfying (7.14). Hence the conclusion follows from Theorem 6.14 extended by Remark 
6.15 — in particular, implication (6.25). □ 

a. Dirichlet case (homogeneous norms). As a consequence of the fact that B^(r2; <f) = 
A^(r2) (Remark 6.12 (i)) we immediately obtain a regularity result in the scale for the 
Dirichlet problem. 

Corollary 7.9. When V = HQ(r2) (Dirichlet problem), assuming that the form a (7.1) is 
coercive on V, if condition (7.14) holds, then any solution of the variational problem (7.4) 
satisfies the regularity result: 

(7.16) fGA^+2(^^) =^ uGA^(n). 

Examplel. 10. For the gradient form ay (for which L is the Laplace operator) on Hj!)(f2),the 
spectrum of the edge Mellin symbol 2le is 

(in 

(l.ll) a(2le) = <^ — , ieZ\{0} 
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and the spectrum of the corner Mellin symbol 2lc is 



(7.18) ^(2l,) = |-l±^;.Di.+ l, 

where /i^^ (for n > 1) is the n-th eigenvalue of the Laplace-Beltrami operator with Dirichlet 
conditions on the spherical cap G^- Hence 



6e(l^,av,Hj(fi)) = - and ay, Hj(fi)) = -1 + X^i; + 1 . 

CJe ^ V ' 4 

b. Neumann case (non-homogeneous norms). For the Neumann problem, it is adequate to 
use the full spaces B^(f2) instead of B^(f2; ^, 0) as in Theorem 7.8, — in the Neumann 
case is empty. These two families of spaces differ by the non-zero Taylor expansions at 
comers for the elements of B^(f2). 

For each comer c the optimal condition on (3^ relating to spaces J™ and B^ relies on the 
condition of injectivity modulo polynomials [19, 16]: The spectrum cr(2tc) has to be replaced 
by the set crj,(2lc) of complex A's for which the condition of injectivity modulo polynomials 
does not hold. This means that there exists a non-polynomial function 

Q 

solution of the problem (7.6) with a polynomial right hand side f on the infinite three- 
dimensional cone /Cc. Note that this condition may differ from the condition in (7.7) only 
for integer A: 

a(2lc)\N = a,(2le)\N. 
Then a, V) is defined as the supremum of the numbers b such that 

{A G C : -i < Re A < 6} n cr^Stc) = 0. 

Theorem 7.11. When V = H^(f2) (Neumann problem), assuming that the form a (7.1) is 
coercive on V, if the following condition holds for the exponents {3 



< -/3e - K h^{^, a, V) Ve G 
-i < -/3c - 1 < &c(^,a,V) VcG<^ 



(7.19) 

' 2 — 2 

then any solution of the variational problem (7.4) satisfies the regularity result: 

(7.20) f G B^+2(fi) ^ u G B^(fi). 

Proof. The proof follows the same steps as for Theorem 7.9, but with one difference: Instead 
of relying directly on [38, Thm 7.1], by a modification of this statement we prove with the 
corner Mellin transform that condition (7.19) implies for any m > max{— /3e, —/3c} the 
regularity u G J™(f2). Once more, the conclusion then follows from Theorem 6.14 extended 
by Remark 6.15 — now, implication (6.26). □ 

Example 7.12. For the gradient form ay on H^(f2), the spectrum of the edge Mellin symbol 
2le is 

(7.21) ^(2i^) = |^![, £gz} 
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and the spectrum of the corner Mellin symbol 2lc is 



(7.22) ^(2ij = |__±^^Neu + _^ ^el^l 

where fi^^^ (for n > 1) is the n-th eigenvalue of the Laplace-Beltrami operator with Neu- 
mann conditions on the spherical cap Gc- Since /i^ i" = 0' the set cr(2lc) contains 0. But one 
can show that the condition of injectivity modulo polynomials is satisfied in A = 0, and that 
it is also satisfied in A = 1 if 1 ^ cr(2lc). Hence we deduce 



6e(fi,av,Hi(fi)) = -^ and ay, Hi(fi)) > min {2, + ^/i.^J + 1} . 

c. A priori estimates along edges. We conclude this section by considerations about the 
nature of necessary and sufficient conditions ensuring the closed range properties along 
edges required by Assumptions 6.5 or 6.13. As mentioned in Remarks 5.5 and 5.13, a 
minimal condition for these assumptions to hold at a chosen edge e G i# is an injectivity 
and closed range condition for the Fourier symbol (l/e(0) ^e,s(07 -^e,s(0) of the system 
{L,Ts,Ds) on the plane sector ICe- 

( Le(0 U = f in /Ce 

(7.23) I re,s(0 u = on Ts, S G ^e, 

[ De,s(0 U = on Ts, S G ^e, 

Here Le is the operator L written in local Cartesian coordinates (x^, xl) G /Ce x IR- The set 
ye is the set of the two faces such that e C s, and for s e S^e, the boundary operators Te^s 
and De,s are the local forms of Tg and D^, respectively. 

• Homogeneous norms. The necessary and sufficient conditions for Assumptions 6.5 to 
hold is that {Le{^) ,Te,siO ^ ^e,siO) defines an operator with trivial kernel and closed range 
from E^^(/Ce) into E°^_,_2(/Ce) for ,^ = ±1. The closed range condition is satisfied if and only 
if,cf[34], 

(7.24) - /3e - 1 ^ Re (a(2le)) := {r/ G M : 3A G (T(2le) with Re A = 7]}. 

In the variational case when the space V is contained in K^^(f2), the trivial kernel condition 
is satisfied as soon as —(3^ — 1 > 0, and even further, for all (3e such that —(3e — 1 > 
—b~{^l, a, V) where &e (fi, a, V) is the supremum of the numbers b > such that 

{A G C : -6 < Re A < 0} n (T(2le) = 0. 

In the Laplace Dirichlet case, the conjunction of the two conditions is 



-/3e-l> and - (3^ - 1 ^ £ E N(£ > 1). 
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• Non-homogeneous norms and Neumann case. Then the necessary and sufficient condi- 
tions for Assumptions 6.13 to hold is that (l/e(0) ^e,s(05 -De,s(0) defines an operator with 
trivial kernel and closed range from J^^(/Ce) into J°^^2(^e) for ^ = ±1. The closed range 
condition is implied by (7.24). The optimal one has to be defined with the injectivity mod- 
ulo polynomials. For the Neumann case, this makes a difference because belongs to the 
spectrum (T(2le) and not to the star spectrum (T*(2le) defined by the injectivity modulo poly- 
nomials. The optimal trivial kernel condition in the Neumann case is —(3^ — 1 > 0. 
More details in the forthcoming work [16, Part III]. 

8. Extensions and generalizations 

In this final section, we describe possible extensions and generalizations of our results. 
More or less straightforward extensions concern non-zero boundary conditions, non-con- 
stant (analytic) coefficients in the two-dimensional polygonal case, and general boundary 
conditions in the polyhedral case. These situations can be handled with the techniques pre- 
sented in the previous sections, and they were omitted here mainly for the sake of brevity. 
Generalizations that could be handled with similar methods, but would need further techni- 
cal work, concern analytic coefficients in the polyhedral case, transmission problems, and 
higher order elliptic systems or systems elliptic in a more general sense. 

8.1. Inhomogeneous boundary conditions and variable coefficients. The fundamental 
estimate of Proposition 2.1 in the smooth case is already formulated in (2.1) for the situation 
of non-homogeneous boundary data. It is also available for variable coefficients. One has 
to introduce the trace spaces on the boundary that correspond to our function spaces on the 
domain. It is well known how to do this, and it is covered in the references given in section 2. 
On the technical side, it is also known how to extend the method of dyadic partitions to the 
case of variable coefficients. Therefore our analytic regularity results in sections 4 and 7 can 
be extended to cover the situation of elliptic systems with analytic coefficients on polygonal 
domains. Such results have been published by Babuska and Guo for a more restricted class 
of elliptic equations, see Remark 7.5. Thus, on a polygonal domain Vt C M^, we can consider 
a general boundary value problem 

{L u = f in f2, 
Ts u = g, on Ts, s G ^ , 
u = hs on r„ s G ^ , 

where L = L(x; d^) is a second order elliptic system, Tg = Ts(x; Sx) and Ds = -Ds(x) are 
boundary operators of order 1 and 0, respectively, and the operators have analytic coeffi- 
cients and may have lower order terms. 

The analogs of Theorems 4.4 and 4.5 are then true, if we augment the regularity assump- 
tions on the right hand side f by the appropriate regularity assumptions on and hg. 

Since Kondrat'ev's results [29] apply to general operators with variable coefficients, the 
basic regularity results for variational solutions are also available, and therefore Theorem 7.1 
can be extended to cover coercive problems with analytic coefficients on polygons. 
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In the three-dimensional case, the extension to non-zero boundary data is as straightfor- 
ward as in the two-dimensional case, but the handling of variable coefficients is of a different 
level of difficulty, due to the anisotropy along the edges. For this, the techniques of analytic 
estimates have to be resumed at a more basic level, involving commutator estimates and 
norms of Sobolev-Morrey type [16, Lemmas 1.6.2 & 2.6.2]. This will be presented in detail 
elsewhere. 

8.2. General boundary conditions. In our theorem on analytic regularity for coercive vari- 
ational problems on polyhedra. Theorem 7.8, we had to restrict the admissible boundary 
conditions to Dirichlet and Neumann conditions. A boundary condition concerning the tan- 
gential components, for example, such as described in (7.3), is not covered by Theorem 7.8, 
although the corresponding result is undoubtedly true. The restriction is not due to the 
tools developed in this paper — the natural analytic regularity shift results in Theorems 6.8 
and 6.14 are proved for solutions of problem (6.3) without this restriction — but due to the 
availability of the basic regularity results that we are quoting from [38], see the proof of 
Theorem 7.8. If one wants to lift this restriction, one therefore has to prove basic regular- 
ity in the appropriate weighted Sobolev spaces for solutions of the boundary value problem 
(6.3). This is outside of the scope of the present paper, but it will be treated in [16]. 

8.3. More general elliptic problems. First we may easily extend the results of this paper 
to transmission problems, namely problem like (6.3) where L has piecewise constant coef- 
ficients (hence some transmission conditions have to be imposed at the common boundary 
of the sub-domains). Indeed an estimate like (2.1) holds for such problems and is proved 
in [16, Theorem 5.2.2]. Second, higher order differential operators like may be treated 
in a similar manner. Finally, our method may be used for the Stokes system (see [25] for 
two-dimensional results). 

The interesting case of boundary value problems for the Maxwell equations is more del- 
icate for several reasons. Whereas the Maxwell equations may be formulated more or less 
equivalently as a second order elliptic system, the boundary conditions will be of a more 
general type than the one treated here. More importantly, the energy space where varia- 
tional solutions are to be found is, in the case of non-convex polygons or polyhedra, not 
contained in H\ and therefore the basic regularity results will be of a different nature. 
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